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ABSTRACT 

We present a comprehensive analysis of structure in the young, embedded 
cluster, NGC 1333 using members identified with Spitzer and 2MASS photometry 
based on their IR-excess emission. In total, 137 members are identified in this 
way, composed of 39 protostars and 98 more evolved pre-main sequence stars 
with disks. Of the latter class, four are transition/debris disk candidates. The 
fraction of exposed pre-main sequence stars with disks is 83%±11%, showing that 
there is a measurable diskless pre-main sequence population. The sources in each 
of the Class I and Class II evolutionary states are shown to have very different 
spatial distributions relative to the distribution of the dense gas in their natal 
cloud. However, the distribution of nearest neighbor spacings among these two 
groups of sources are found to be quite similar, with a strong peak at spacings 
of 0.045 pc. Radial and azimuthal density profiles and surface density maps 
computed from the identified YSOs show that NGC 1333 is elongated and not 
strongly centrally concentrated, confirming previous claims in the literature. We 
interpret these new results as signs of a low velocity dispersion, extremely young 
cluster that is not in virial equilibrium. 
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Introduction 



Observations of embedded, star-forming clu sters and groups show that the stellar distri- 



butions are often elongated, clumpy, o r both (cf. iLada fc Lada 



2003; Gutermuth et al. 2005 



Teixeira et al.ll2006l : lAUen et al.l 120071 ). and the structure seems tied to the distribution of 
dense gas in the clusters' natal molecular clouds. Most clouds have some non-spherical 
structure, yet current results suggest that the associated clusters have varying degrees of 
agre ement with the cloud's structure depending on how deeply the members are embed- 
ded (IGutermuth et al.ll2005l : iTeixeira et al.l 120061 ). Given the high frequency of asymmetric 
structure in clouds, it seems reasonable to assume that the exposed, relatively structureless 
clusters we observe may have been asymmetric in a previous epoch. With the ejection of the 
majority of their natal gas and adequate time to migrate from their birth sites, the imprint 
of the underlying cloud structure could be lost rather quickly in a recently exposed young 
cluster. Thus the structure we measure in a distribution of young stellar objects (YSO) 
relative to the dense gas in the associated cloud may be a reasonable proxy for the current 
dynamical state of an embedded cluster. 

The launch of Spitzer has provided a potent new facility for studies of the structure 
of young and embedded clusters. While some members of these young clusters are disk- 
less pre-main sequence stars (Class III) , the majority of t he membership are sources with 
excess emission at mid-IR wavelengths (IHaisch et al.ll200ll ). made up of a mixture of proto- 
stars (Class I), still embedded and accreting from dense spherical envelopes, and the slightly 
more evolved pre-main sequence stars with circumstellar disks (Class II). The mid-IR spec- 
tral energy distributions (SED) of these objects are dominated by the emission from their 
dusty circumstellar material, making them easily distinguishable from pure photospheric 
sources such as unrelated field stars and indistinguishable diskless cluster members. Spitzer's 
sensitivity to mid-IR emission makes it the best tool currently available for identifying 
and characterizing YSOs with IR-excess, and that sensitivity is sufficient to detect these 
sources down to the Hy drogen-burning mass limit for regions within the nearest kiloparsec 
( IGutermuth et al.l 120041 ). The complete census of YSOs with disks in a young cluster repre- 
sents a high-confidence sample of bona fide cluster members, and for t he youngest regions , 
such a sample includes a high fraction of the total number of members (IHaisch et al.ll200ll ). 
Furthermore, by separating the two canonical evolutionary classes of YSOs that Spitzer so 
effectively detects, we are able to probe both the recent overall star-forming activity in the 
region as traced by the stars with disks, and the immediate star formation as represented 
by the protostars, since this phase is expected to be short-lived relative to the former. 



Many young clusters are dominated by heavy and spatially variable extinction from dust 
within their natal molecular cloud environment. Another advantage Spitzer hhngs to studies 
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of these regions is that the mid-IR wavelengths targeted by Spitzer's imaging instruments 
(Infrared Array Camera, or IRAC, at 3.6-8.0 /im and Muhiband Imaging Photometer for 
Spitzer, or MIPS, at 24-160 fim) are l ess affected by extinction from du st in comparison 



to near-IR or visible wavelengths (e.g. ICardelli. Clayton, fc Mathid Il989l ). Recent papers 



have done an excellent lob of characterizing the reddening law in the IRAC bandpasses 



(llndebetouw et al.l 12005 



Flaherty et al.l 120071 : iHuard et al. 2007 ) and have made a first at 



tempt at doing the same for the 24 /im channel of MIPS ( iFlaherty et al.l 120071 : iHuard et al. 



20071 ). Given the wealth of information Spitzer can bring to the study of embedded clusters, 
we have surveyed over thirty clustered star-forming regions within the nearest kiloparsec 
through the Guaranteed Time Observations (GTO) program for the IRAC and MIPS in- 
strument teams. With these data, we can not only provide a full census of the YSOs with 
IR-excess in each region, but also perform a detailed examination and comparative analysis 
of structure in young, embedded clusters. 

NGC 1333 has been a popular target for observations of deeply embedded protostars 
via radio wavelengths due t o the presence of an unprecendented number of molecular out- 
flows (e.g. iKnee fc Sandelll 120001) assoc iated with several bright IRAS sources considered 
Class protostars (I Jennings et al.lll987l ). all in relatively close proximity to the Sun (250 pc, 
Enoch et al.ll2006l ). Many of the outflows are traced by shock- induced em ission, a clear sign 



that they are indeed affecting the local, quiescent cloud material (e.g. IWalawender et al. 



20051 ) . Because of this, some studies have claimed that the NGC 1333 den se molecular cloud 



core is in the process of being destroyed b y influence from outflows (e.g. iWarin et al.lll996 
Sandell &: Kned l200ll : iQuillen et al.ll2005l ). In addition to the numerous protostars, there 
is a cluster of pre-ma. i n sequence stars identified by number counts analysis in the near - 



IR (Strom et al 



1976 



Aspin et all ll994J : iLada. Alves. fc Ladal Il996l : IWilking et al.l l2004h . 



Lada. Alves. fc Lada ( 19961 ) suggested that the distribution was well-described as a "double 
cluster" , having two distinct surface density maxima. 

Here we present a Spitzer IRAC and MIPS imaging and photometric analysis of the 
NGC 1333 young cluster, one of the most nearby large membership (A^ > 100) clusters 
in the Spitzer Young Cluster Survey. We achieve a near-complete census of the cluster 
membership that possesses circu mstellar material, si gnificantly surpassing the sensitivity of 
ground-based mid-IR surveys (cf. iRebuU et al.ll2003l . and references therein.). In addition, 
we statistically infer the population of pre-main sequence stars that lack disks, enabling 
an estimate of the fraction of members with disks. Finally, using the identified YSOs, we 
apply both established and recently developed methods for characterizing the structure of 
this cluster, in specific reference to previous claims of structure found in the literature. 
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Observations 



2.1. Multi-band Point Source Photometry 



2.1.1. IRAC Imaging 



Observations of th e NGC 1333 region y ere obtained with the IRAC instrument (IFazio et al 



20041 ) onboard Spitzer flWerner et al.ll2004l ) over two epochs: February 10, 2004 (PID 6) and 
September 8, 2004 (PID 178). The total area covered by all four IRAC bandpasses in the final 
mosaic is roughly 40' x 30', centered on RA(J2000) = 03^29'"00^ Dec.(J2000) = 31'^20™00^ 
The typical total exposure time per band is 41.6 seconds, comprised of four 10.4 second 
integrations. High Dynamic Range (HDR) mode was used, so there are matching images 
at 0.4 second integration times to obtain photometry for bright sources. Mosaics were built 
at the native instrument resolution of 1.2" per pixel with the standard Basic Calibrated 
Data (BCD) products provided by the Spitzer Science Center from their standard data pro- 
cessing pipeline version S13.2. Additional processing of bright source artifacts ("jailbar", 
"pulldown", "muxbleed", and "banding") is performed on these data via customized ver- 
sions of Interactive Data Language ( IDL) scripts developed by the IRAC instrument team 
(IHora et al.l |2004J : iPipher et al.l 120041 ). Mosaicking was performed using Gutermuth's WC- 
Smosaic IDL package, which integrates the following features: a redundancy-based cosmic 
ray detection and rejection algorithm; frame-by-frame distortion correction, derotation, and 
subpixel offsetting in a single transformation (to minimize data smoothing); and on-the- 
fiy background matching. WCSmosaic is built with heavy dependence on the FITS and 
WCS access and manipulation procedures provided by the IDL Astronomy Users Library 
(jLandsmarull993l ). 



Point source detection and aperture photometry is performe d via Gutermuth's IDL pho- 
tometry and visualization tool. Phot Vis (iGutermuth et al.ll2004l ). version 1.10. Of significant 
note in this new version of Phot Vis is the inclusion of a spatially variable source detection 
threshold made possible through the computation of a median absolute deviation uncer- 
tainty map from the Gau ssian-convolved image computed as part of the DAOfind detection 
algorithm (jStetsonI 119871 ). With this modification, the improved algorithm is able to more 
clearly distinguish point sources from background variations, dramatically improving source 
detection robustness in regions of bright, structured nebulosity common in IRAC images of 
active star-forming regions. The full-width at half-maximum (FWHM) of every detection 
is measured and all detections with FWHM > 3.6" are considered resolved and removed. 
Some marginally resolved but highly structured shock emission sources were also visually 
identified and removed. Aperture photometry was performed on all sources with an aper- 
ture radius of 2.4" and background flux estimation was done within concentric sky annuli of 
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i nner and outer radi i of 2.4" and 7.2" respe ctively, following the method of previous work 
( iMegeath et al.ll2004j : iGutermuth et al.ll2004l). Photorn etric zero points were derived directly 
from the calibrations presented in iReach et al. ( 2005 ). using standard aperture corrections 
for the radii adopted. See Table [T] for 90% differential completeness limits and zero-point 
magnitude calibration values adopted for all bandpasses and all data sources. Phot Vis is 
built with heavy dependence on the FITS access and manipulation pro cedures and DAO phot 
version 1 port to IDL provided by the IDL Astronomy Users Library (IL andsmani 1 1 9 93l ) . 



2.1.2. MIPS Imaging 



MIPS flRieke et al.ll2004[ ) imaging of NGC 1333 was obtained on February 3, 2004 (PID 
58) at the medium scan rate. Twelve scans of 1.75 degrees at 160" offsets were needed to 
cover the full IRAC coverage. Image treatment and mosaicking was performed with the 
MIPS instrument team's Data Analysis Tool (DAT). The DAT is used to treat bright source 
artifacts, persistence artifacts, and cosmic ray hits in each image. It then mosaics the images 
together, removing spatial scale distortion before integrating each frame into the mosaic at 
the native MIPS resolution of 2.54" per pixel. Due to the low resolution and sensitivity of 
the 70 /xm and 160 /im data obtained simultaneously, only the 24 /im data are presented 
and used in this paper. An ana lysis of these data for the entire Pe rseus molecular cl oud is 
presented in iRebuU et al.l (120071 ) as part of the c2d Legacy Survey (lEvans et al.ll2003l ). 



MIPS 24 /im point source detection and aperture photometry were also performed with 
Phot Vis version 1.10. Aperture and inner and outer sky annulus radii were chosen to be 
7.6", 7.6", and 17.8", respectively. The photometric zero point of [24] = 14.6 for a 1 DN/s 
flux measurement with the adopted aperture and sky annulus radii was chosen to make the 
mean [5.8] — [24] colors of the nine photospheric sources detected with MIPS equal to zero 
(reduced from thirteen sources after applying iterative statistics techniques). RMS scatter 
of the internal measurement is 0.04 magnitude. 



2.1.3. 2MASS JHKs Photometry 

All sources detected by Spitzer are matched to Two-Micron All Sky Survey (2MASS) 
Point Source Catalog (PSC) sources, extending the potential wavelength coverage for these 
objects to J band at 1.2 fim. The NGC 1333 Spitzer survey region is also one of the 
locations covered by deeper 2MASS observations, and the merged data from these deeper 
observations were kindly provided by Roc Cutri ahead of public release. The deeper data 
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is reduced usin g the 2MASS final release data pipeline, which is documented in detail in 



Skrutskie et al.l (120061 ) . These data are taken from observations of six times the integration 



time of the standard 2MASS observations as part of the 2MASS extended mission, yielding 
approximately a 1 magnitude improvement in the data completeness limit and reducing 
measurement uncertainty for sources within a magnitude of the detection limit of the regular 
2MASS data JCarpenter et al.lboOlh . 



2.1.4- The Final Catalog 

Source list matching for the final catalog is performed in stages. First, the four IRAC 
band source lists are matched using Gutermuth's WCSphotmatch utility, with a maximum 
1" radial tolerance among matches. Second, the 2MASS catalog is integrated with the same 
software, using the mean positions of all entries in the merged IRAC catalog. The RMS 
scatter of the radial residual offsets between each IRAC band and 2MASS ranges from 330 
to 390 milliarcseconds. Finally, the MIPS 24 fim catalog is integrated with WCSphotmatch 
using a wider 3" maximum radial tolerance, now using the mean 2MASS and IRAC positions 
for all entries in the second stage catalog. The RMS scatter of the radial residual offsets 
between the MIPS catalog and 2MASS is 750 milliarcseconds. 



2.2. Molecular Cloud Maps 



Tracers of the distribution of gas in a molecular cloud have varying levels of reliability, 
depending on environment, and that of NGC 1333 is especially challenging. Many recent 
papers have shown that CO and its isotopes are relatively poor column density tracers (e.g. 
Kramer et al.l Il999l ) outside of a rather narrow range of column densities (A. Goodman, 
priv. comm.), and the many active outflows in NGC 1333 add considerable confusion. The 
high galactic latitude of NGC 1333 has low field star surface density and a general lack of 
signi ficant numbers of background giants compared to regions like the galactic bulge (e.g. 
B68, lAlves et al.ll200ll ). thus extinction maps derived from the reddened near-IR colors of 
background sources are generally low resolution. In addition, the high overall extinction 
{Ay > 10) in NGC 1333 and the presence of a partially embedded population of pre-main 
sequence stars will bias the map toward low er column densities, pa rticularly where cluster 
members are more numerous than field stars (iGutermuth et al.ll2005l ). Both cold (e.g. in the 
cold dust filament where the Class objects IRAS 4a, 4b, and 4c are forming) and hot (e.g. 
in the northern reflection nebulosity illuminated by late B stars) environments are found 
in NGC 1333, making conversion of submillimeter dust continuum emission to gas column 
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density particularly uncertain (jSchnee et al.ll2005l ). Recent work y y ith ot her dense gas tracers 
holds promise, such as the N2H~^ map presented in IWalsh et al.l (120061 ) . but this particular 
study (11' X 11') does not cover the large field of view of the Spitzer imaging data. 

Given the many pitfalls, we adopt two different cloud mapping methods here based only 
on their utility for two specific functions: tracing large scale, low-level extinction reliably, and 
tracing dense gas morphology reliably. For the former purpose, we create extinction maps 
using the H — Ks colors of field stars from th e combined 2MASS data . The method used 
is identical to that presented in Sec tion 3.1 of iGutermuth et al. I J2OO5I), based on elements 
of both the NICE Jhada et al.lll994[ ) and NICER Jlombardi fc AlveslboOlh algorithms. The 
resulting maps are sensitive to extinctions from 1 Ay to 18 Ay and have variable resolution 
based on local stellar surface densities, ranging from 70" to 200". By characterizing low- 
level extinction well, we are able to customize our color- and magnitude-based cuts to 
mitigate extragalactic contamination to utilize the fact that extinction shifts the number 
density of contaminators per magnitude to dimmer magnitudes (see Section WA\ below) . For 
tracing the morphology of dense gas, we adopt the JCMT/SCUBA 850 /im map built from 
the combination of all such observations of this region found in the JCMT online archive 



(iKirk et al.l 120071). Th e bulk of these observations were originally pub lished and analyzed i n 



Sandell fc Kneel (12001 ) with expansion of the original field presented in Hatchell et al. ( 2005 ). 
While derived column densities from this map are highly uncertain, 850 /xm fiux should trace 
the overall morphology of the densest portion of the cloud with reasonable accuracy to enable 
a comparison with the spatial distribution of cluster members (see Section EH] below) . 



3. The Spitzer View of NGC 1333 

We present color renditions of the Spitzer IRAC and MIPS mosaics in Figures [1] & [^. 
In Figured], we have mapped the IRAC 3.6, 4.5, and 8.0 //m mosaics in blue, green, and red 
respectively. NGC 1333 is unique in the Young Cluster Survey in that it is dominated by a 
complex network of 4.5 /xm-dominated (dark green) shock emission from several different out- 
flows colliding with ambient cloud material. A more diffuse, blue-green color shows regions 
of scattered light from the surfaces of the molecular material. Finally, the bright red region 
in the north coincides with the visibly exposed reflection nebulosity that is "NGC 1333". 
The 8.0 /zm emission in red here is commonly considered to be UV-excited polycyclic aro- 



For reference, the near-IR view presented in Fig. 16 of lLada k Ladal (|2003h is the deepest near-IR image 
of the entir e main cluster regio n currently available, and many of the region's well-studied objects are marked 
in Fig. 1 of IWalsh et all |2006l) . 
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matic hydrocarbon (PAH) feature emission, where the UV source is usually one or more 
intermediate or high-mass stars. There is also a less intense, but similarly sized region of 
PAH emission in the south end of the region, with the IRAS 4abc molecular cloud filament 
seen in absorption (RA = 03:29:07, Dec. = +31:13:30). This is a new feature revealed by the 
extinction-penetrating capabilities and mid-IR sensitivity of Spitzer, likely a bubble being 
evacuated by a previously unknown massive star. We have identified a candidate for the 
exciting star, ASR 54, located near the projected center of the newly revealed PAH neb- 
ulosity structure and apparently behind the dense, dusty molecular cloud filament. While 
not detected in th e 2MASS dataset presented h ere, ASR 54 was dete cted at i^' = 15.74 by 



Aspin et al.l (Il994j ). Given the extinction law of iFlaherty et al.l (120071 ). the IRAC photome- 



try of ASR 54 is consistent with a heavily extinguished photosphere [Ak ~ 10). As such, for 
an assumed distance of 250 pc, we estimate that the absolute i^-band magnitude of ASR 54 
is K = —1.25, consistent with a late B-type star. The 24 /im photometry of ASR 54 also 
suggests there is some excess emission in this bandpass, resulting in this source's inclusion 
as a transition/debris disk candidate (see below). 

In Figure Ej we have mapped the IRAC 4.5, 8.0, and MIPS 24 /xm mosaics in blue, 
green, and red respectively. In these colors, the many protostars in this region are made 
visually identifiable by the bright red halo around each source. The pre-main sequence stars 
with disks and diskless stars are noticeably bluer in this view. The bright red source on the 
east side of the cluster in this image (RA = 03:29:37.1, Dec. = +31:15:46) is detected only 
at 24 yum. Given the elongated nature of the detection, no detections in any other IRAC 
or MIPS bands, and its position offset from the molecular cloud core, we are confident that 
this source is an asteroid. 



4. The Census of YSOs in NGC 1333 

Here we describe the empirically-derived methods used to classify sources with infrared 
excess in NGC 1333. We have adopted a three phase approach to achieve maximum sensitiv- 
ity to YSOs in all evolutionary states that exhibit IR-excess under a wide range of nebulosity 
environments from the combined multiband photometric dataset. In each phase indepen- 
dently, we account for and attempt to mitigate contamination from extragalactic sources 
that masquerade as bona fide YSOs in the field. 
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4.1. IRAC Four-Band Source Characterization 



As our primary dataset, we consider tlie set of all sources with detections (cr < 0.2 mag) 
in all four IRAC bandpasses, a total of 578 sources. There are now several methods in 



the literature to identify and classif y YSOs with such photometry (cf. lAUen et al.l 12004 



Megeath et al.|[2004 : Lada et aPboOG ). Each of the cited schemes depends in large part on the 
3.6 fim photometry to distinguish between YSOs of differing evolutionary class. However, in 
deeply embedded clusters such as NGC 1333, YSO classification schemes that rely on 3.6 /im 
photometry will suffer to some degree because of the degeneracy caused by reddening from 
dust. It has been shown that the extinctio n law flatte ns considerably through the 4.5-8.0 fim 
regime, only rising in t he 3.6 /xm bandpas s ( Lutzll999h. and this has been confirmed in se veral 



Spitzer-hased studies ( llndebetouw et al.l 12005 



Flahertv et al.ll2007l : iHuard et all 120071 ). To 



dramatically reduce the risk that differential extinction will confuse our classification of 
YSOs, we have adopted a new scheme that relies primarily on [4.5] — [5.8] color to distinguish 
between Class I and Class II YSOs. This color is much less affected by dust extinction than 
any [3.6]-based color , yet it can still be emp loyed to distinguish between Class I and Class II 
YSOs (see Fig. 2 of iHartmann et al.ll2005[ l. 



Before attempting to identify YSOs from these data, we must first remove various 
extragalactic contaminants that may be misidentified as YSOs. Initially, we identify and 
remove 13 sources dominated by PAH feature emission, likely star-forming galaxies and 
weak-line active galactic nuclei (AGN), following the method described in Appendix |A] (see 
the color-color diagrams in Figure [3]). We then identify and remove 31 sources that fall 
within the region of color-magnitude space dominated by broad-line AGN (see Figure H]). It 
is important to note here that the presence of a molecular cloud introduces a bias against 
a source of a given magnitude being an AGN, as extinction from dust in the cloud shifts 
the density of extragalactic sources per magnitude toward dimmer magnitudes. To account 
for this, we artificially assume that every source is located behind the cloud and deredden 
their photometry using the value from our extinction map at each source's position. 
Then we apply the AGN cut in [4.5] vs. [4.5] — [8.0] color-magnitude space described in 
Appendix [A] to the dereddened photometry. Those sources flagged in this way have colors 
entirely consistent with low luminosity YSOs, thus we will reexamine these sources for signs 
of a protostellar signature using the MIPS 24 fim photometry in Section H73] below. See the 
left plot of Figure O for the spatial distribution of sources filtered from the NGC 1333 source 
list as likely contaminants. 

To test the filtration process, we compare the resulting distribution of PAH emission 
and AGN candidate sources with models derived from our re-reduction of the Spitzer survey 
for galaxies in the Bootes field (see Appendix |A]) . These models are built from the [5.8] 
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magnitude distribution shown in Appendix \^ scaled for the area of the NGC 1333 IRAC 
four-channel coverage and adjusted to account for the relative areas of differing extinc t ion vi a 



the 2MASS-derived extinction map and the IRAC reddening law of iFlaherty et al.l (120071 ) . 
See Figure M for comparison histograms for the PAH and AGN sources. Overall, the agree- 
ment between the filtered source distribution and the models is very good, suggesting that 
very few bona fide YSOs were mistaken as contaminants and discarded. However, there is 
likely residual contamination of 1 ± 1 sources after application of the contaminant filtration 
process. According to the Bootes characterization, there is a higher probability that these 
are mistaken as Class I rather than Class II sources (see Appendix E|) . 

We now attempt to filter one final contaminant from the source list before we proceed 
with identifying and classifying the YSO population: unresolved blobs of shocked emission 
from high-velocity outflows interacting with the cold molecular cloud. These sources have a 
particularly large 4 .5 /im excess due to shocked molecular hydrogen line emission in this band 



( ISmith et al.l 120061 ). and this fact is exploited to pull these sources out of the final list. As 
was the case for the list of likely AGN, the shock emission source list can hold some bona fide 
protostars. Hence we also will check this list for MIPS detections in Section IT3l NGC 1333 
has the most structured shock emission of any cluster in the Spitzer Young Cluster Survey, 
and thus we have used it to empirically determine the region of IRAC color space that these 
sources occupy. Based on our findings, all sources with photometry that obeys all of the 
following constraints are likely dominated by shock emission and thus are removed: 

[3.6] - [4.5] > ^ X (([4.5] - [5.8]) - 0.3) + 0.8 

[4.5] - [5.8] < 0.85 
[3.6] - [4.5] > 1.05 

Four sources are removed using these contraints. The remaining sources are dominated 
by stellar sources, and infrared excess sources still present are true YSOs of both Class II 
and Class I evolutionary states. As mentioned above, we have adopted the [4.5] — [5.8] 
color as the primary discriminant between protostars and stars with disks in this empirically 
determined classification scheme. It has been shown in the low extinction environment of 
the nearby Taurus mol ecular cloud that YSOs can be successfully classified with this color 



(IHartmann et al.ll2005l ). and the reddening effect of dust along the li ne of sight is much 



less severe for this color than for [3.6] — [4.5] (e.g. iFlaherty et al.l 120071 ). Sources are likely 
protostars if they have extremely red discriminant color ([4.5] — [5.8] > 1). Additionally, 
any sources with moderately red discriminant color (0.7 < [4.5] — [5.8] < 1.0) that also have 
[3.6] — [4.5] > 0.7 are likely protostars, though in rare cases, a highly reddened Class II source 
could have these colors as well. 



Once the protostars are removed from consideration, we move on to identifying the 
Class II sources from those that remain. All sources with the following color constraints are 
likely Class II sources: 



[4.5] - [8.0] > 0.5 



[3.6] - [5.8] > 0.35 



[3.6] - [5.8] < ^ X (([4.5] - 



8.0]) -0.5) + 0.5 



In total, we identify 29 sources with protostellar colors and 93 with colors consistent 
with a pre-main sequence star with a circumstellar disk using this method. All other sources 
still remaining are consistent with photospheric colors. Almost all of these objects are field 
stars unrelated to the cluster itself. Unfortunately, pre-main sequence stars that lack disks 
but are bona fide members of the cluster are still mixed in with this group and cannot be 
identified from the data presented here. For this work, we simply account for these sources 



statist ically, using star counting methods similar to those presented in iGutermuth et al. 
fl2005h (See Section [6] below) . 



This classification scheme has been verified using the JHKs and I RAC photometry 



for th e nearby, low-extinction Taurus molecular cloud YSOs presented in iHartmann et al. 
(120051 ). Of the 82 Taurus sources presented in that work, 72 have full four band IRAC 



photometric coverage with a < 0.2 mag. Our method classifies as Class III all 24 source s 
considered Class III in that work. Of the 41 classified as Class II in iHartmann et al.l (120051 ). 
our method classifies 39 as Class II and 2 as Clas s III. Finally, we c lassify as Class I 4 of the 
6 sources called Class 0, Class I, or Class I/II in IHartmann et al.l (120051 ): the other two we 
classify as Class II. 

The discrepant sources in each evolutionary class have measurably weaker IR- excess than 



the no n-discrepant ones, as measured by the aiRAc classification described in iLada et al. 
(120061 ). In that work, the slope of the linear least squares fit to the 3.6-8.0 /xm SEDs in 



log \F\ vs log A space, defined as aiRAC^ was used to separate sources by evolutionary 
class. In t he case of the two sou rces we have classified as Class III but are classified as 
Class II in IHartmann et al.l (120051 ) . t he ajRAc v a lues a re -3.07 and -2.66, both consistent 



with photospheric SEDs according to iLada et al.l (l2006l) . Similarly, the two sources we call 
Class II that are considered Class I in IHartmann et al.l (120051) have amAc values of -1.59 



and -1.32, consistent with Class II sources (ILada et al 



20061 ) . Regardless of the ambiguity 



of the classifications for these few sources, this comparison demonstrates the overall efficacy 
of the new classification technique. 
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4.2. Additional YSOs Identified via a JHKs[S.6][4:.5] YSO Classification Scheme 

NGC 1333 is a moderately dense cluster with significant structured nebulosity (scat- 
tered light, H2 shocks, and PAH emission), and both high stellar density and bright, diffuse 
emission can reduce local source sensitivity. This issue is made worse by the decreasing 
angular resolution, large and bright PSF wings, and intrinsically lower sensitivity of the 
longer wavelength mosaics from the 5.8 and 8.0 /im IRAC channels and MIPS 24 fim. As 



mentioned above, i t has been shown that [3.6] 



discriminant color (lAUen et al.l |2004J : iMegeath et a. 



.5] is a useful YSO evolutio nary class 



2004 



Hartrnann et al.ll2005l ). but suf- 



fers from degeneracy with differential extinction. iGutermuth I (120051 ) utilized dereddened 
2-8 /im SED slopes to demonstrate that one could separate protostars, stars with disks, 
and diskless photospheres quite effectively with the K — [3.6] vs. [3.6] — [ 4.5] color-color 



diagra m, provided that the colors were dereddened via the reddening law of [Flaherty et al. 
J2007h . This reddening law was adopted because it was consistently measured through sev- 



eral nearby, star- forming molecular clo uds, and shown to be si gnificantly different than the 
reddening law from diff use ISM dust (llndebetouw et al.ll2005l ). This empirical method of 
YSO identification from Gutermuth ( 20051 ) is explained below. 



First we measure the line of sight extinction to each source as represented by the ^ 



co lor excess rat i o, uti lizing baseline colors based on th e Classical T Tauri St ar (CTTS) locus 
of iMeyer et al.l (119971 ) and standard dwarf star colors (IBessell fc Brettlll988l ). To accomplish 
the latter, we force [J — H]q > 0.6, a simplifying approximation for the intrinsic colors of 
low mass dwarfs. Here are the equations used to derive the adopted intrinsic colors from the 
photometry we have measured: 



[J - H]o = 0.58 x[H- K]o + 0.52 



[H - K]o = [H-K] 
[H - K]o 



i[J-H] 



[J-H]o) X 



H-K 



Ej-H 

[J ~ -f^]meas ~ Eh~^ ^ ~~ K\meas ~ 0.52 

0.58 



Eh- 



Once we have measured the component of the H — K color excess that is caused by 
reddening, we compute t he dereddened K — [3 .6] and [3.6] — [4.5] colors using the color 
excess ratios presented in iFlahertv et al.l (I2007D . specifically = 1.73, = 1.49, 



and 



Eh- 



X-[3.6] 



i'K-[4.5] 



1.17: 



[K - [3.6]]o = [K - [3.6]], 



{[H-K]meas-[H-K], 



E 



X 



K-[Z.6\ 



E 



H~K 
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[[3.6] - [4.5]]o = [[3.6] - [4.5]] 



meas 



{[H-K] 



meas 



[H - K]o) X 



Eh-k 




n -1 




-1 



-5'[3.6]-[4.5] 

Eh-k 



We now identify additional YSOs as those sources with dereddened colors that obey all 
of the following constraints: 



Among these YSOs, sources that also follow the following selection are considered likely 
protostars (the rest are considered Class II): 



To minimize the inclusion of dim extragalactic contaminants, we apply a simple bright- 
ness limit using the dereddened 3.6 yum photometry. All sources classified as Class II with 
this method must have [3.6]o < 14.5 and all protostars must have [3.6]o < 15. In total, the 
use of this supplemental method allowed us to consider 504 more sources, and we added one 
Class I and four Class II sources to the NGC 1333 YSO list. 



4.3. Adding and Checking YSOs with MIPS 24 //m Photometry 

We utilize 24 ^m. photometry from the MIPS instrument in several ways. First, sources 
that were considered photospheric in the previous two phases but have excess emission at 
24 yum ([5.8] — [24] > 1.5) are very interesting. These sources are a mixture of regenerative 
debris disks and "Transition Disks"; the latter are Class II sources with significant dust 
clearing within their inner disks. We have required that valid sources in this group also have 
[3.6] < 14 in an effort to minimize contamination from extragalactic contaminants. 

Second, the NGC 1333 cluster is well known partly because there arc many known 
Class protostcllar members driving spectacular outfiows. Many of these sources are so 
deeply embedded that we cannot extract usable photometry in one or more of the four 



[[3.6] - [4.5]]o > 0.101 
[K - [3.6]]o > 
[K - [3.6]]o > -2.85714 x ([[3.6] - [4.5]]o - 0.101) + 0.5 



[K - [3.6]]o > -2.85714 x ([[3.6] - [4.5]]o - 0.401) + 1.7 
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IRAC bands. As such, we note any source as a likely deeply embedded protostaii if it 
lacks detections in some IRAC bands yet has bright MIPS 24 /xm photometry ([24] < 7 
and [X] — [24] > 4.5 mags, where [X] is the photometry for any IRAC detection that we 
possess). The adopted magnitude limit at 24 /im is needed to mitigate extragalactic source 
contamination. The specific value th at we adopt here is slightly more conservative than that 



adopted by c2d (IHarvey et al.ll2006l ). With the addition of the 24 /xm data, we have more 
information to filter lower luminosity YSOs from our AGN candidates and shock emission 
dominated sources flagged in Section 14. 1[ We reinclude flagged sources as likely protostars 
if they have both bright MIPS 24 fim photometry ([24] < 7, as before) and convincingly red 
IRAC/MIPS colors ([3.6] - [5.8] > 0.5 and [4.5] - [24] > 4.5 and [8.0] - [24] > 4). 

Finally, all previously identified protostars that have 24 /im detections are checked to 
ensure that their SEDs do indeed continue to rise from IRAC to MIPS wavelengths. All 
protostars that have MIPS detections must have [5.8] — [24] > 4 if they possess 5.8 fim 
photometry, otherwise they must have [4.5] — [24] > 4. If a source does not meet this color 
requirement, it is likely a highly reddened Class II. Of the 29 sources identified as protostars 
above, we find only one source that was mistakenly classified as a protostar. We reclassify 
this source as a Class II accordingly. 

In total, we identify four transition/debris disks using the above integration of MIPS 
24 /im data. Of the 99 sources with MIPS detections and incomplete IRAC photometry, 
we add 11 deeply embedded sources to our list of protostars. Many of the latter deeply 
embedded sources are known protostars from the literature, such as IRAS 4a, 4b, and 4c. We 
also reclassified one deeply reddened Class II that was first classified as Class I in Section H]T] 
above. 



4.4. YSO Census Summary 



For the entire three phase method, we find 4 transition/debris disks, 94 Class II sources, 
and 39 protostars, for a total membership with disks of 137 sources in NGC 1333. For 
comparison to other studies, see the top plot of Figure [8] for the [3.6] — [4.5] vs. [5.8] — [8.0] 
diagr am used in sever a l recent Spitzer-hased papers to identify an d/or characterize YSOs 



le.g- 



Allen et all |2004J : iMegeath et all |2004J : iHartmann et al.ll2005l ). The photometry and 



classifications are presented in Tables [2] & [3] (Note that the entries and indexing of these 



^Classifying sources as bona fide Class protostars requires a detailed treatment including submillimeter 
flux measurements which are beyond the scope of this paper. As such, these sources are still labelled Class 
I in Table [3l but with the addition of an asterisk to note their anomalous nature. 
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tables are in ascending Right Ascension order /it separated by YSO class as documented 
in Table O). For the rest of this paper, we consider the transition/debris disks as part of 
the Class II group, yielding a total of 98 Class II sources for the analyses that follow. The 
spatial distribution of these sources is shown in the right plot of Figure O While the sources 
are indeed quite concentrated compared to the extragalactic contaminants in the left plot of 
Fig. [5], the structure of the NGC 1333 cluster is visibly non-uniform. 



5. Spatial Distributions of YSOs by Class 



With a confident census of the YSOs in this region, we now examine the spatial distribu- 
tion of these sources in relation to each other and that of the dense natal gas. In Figure [TOb. 
we have plotted the distribution of known pre-main sequence stars, namely the Class II and 
transition/debris disk sources, on contours of 850 fim dust continuum emission. This popu- 
l ation of sources reflects the "double cluster" distribution reported by iLada. Alves. fc Lada 
(119961 ). The two main local peaks in stellar density are confined to local minima in the 
dust emission. In contrast, we have plotted the distribution of likely protostars on the same 
contours in Figure [TOb. Not surprisingly, these sources trace the dense gas more closely, and 
thus they bear little resemblance to the "double cluster" distribution. 

Overall, the two-dimensional morphology of the combined star-gas system is strongly 
filamentary, and bears the appearance of a central elongated ring with three filaments ex- 
tending away from the ring radially, roughly evenly spaced in azimuth. The nodes where 
these filaments meet the central ring structure seem occupied by many young stars, while 
the center of the ring has a noticable absence of stars and gas. A ring -like morphology in 
dense gas has been detected before in the Coalsack (ILada et al.l 12004), and some rnodels 



of collapsing clouds have y ielded various ring-like morphologies (e.g. iLi &: Nakamural 12002 
Burkert fc Hartmanrul2004l ). Thus, it is possible that the appearance of ring-like structure 
in NGC 1333 is not a spurious effect of multiple filamentary structures partly overlapping 
along the same line of sight, but this possibility cannot be ruled out with these observations. 



5.1. Source Spacing and Proximity to Dense Gas 

Here, we analyze the spatial distributions of the separate YSO populations and how they 
relate to the dense gas in a quantitative way. Figure [11] is a histogram of the SCUBA 850 /xm 
fluxes at each source's position for the Class I and Class II sources separately. As demon- 
strated qualitatively in Figure [101 protostellar sources are much more strongly confined to 
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the regions of dense gas traced by the submiUimeter emission than their more evolved Class II 
counterparts. From these data, we choose a flux threshold of 0.05 Jy/beam to demonstrate 
that most (27/39 = 69.2%) Class I sources trace dense gas, while a considerably smaller 
fraction of Class II sources (20/98 = 20.4%) are found in regions of such high submiUimeter 
flux. 

However, while their relationship to the dense gas is quite different. Figure [12] shows 
that the nearest neighbor spacings between sources are quite similar for each evolutionary 
class, as well as for the entire set of YSOs combined. There is a clear peak in each histogram 
between 0.024 pc (5000 AU) and 0.048 pc (10000 AU). A similar peak has recently been 
found in the nearest neighbor spacings of the "Spokes" cluster o f protostars in NGC 2264, 



where most protostars are separated by 12000 AU to 20000 AU (ITeixeira et al.ll2006l ). The 
peak in nearest neighbor spacings among the NGC 1333 protostars and stars with disks is 
half of the value found for the "Spokes" cluster. However, this peak is consist ent with the 



clump and YSO separations of ~6000 AU observed in the Ophiuchus region by lMotte et al. 



(]l998h . 



5.2. Surface Density Mapping 



Following the method presented in lGutermuth et al.l (120051 ) . we have constructed surface 
density maps for the identified YSO distribution in NGC 1333 using the Nearest Neighbor 
methocEI. Figure [T3l shows a comparison between maps generated using n = 4, 6, 11, and 18 
neighbors. The smallest n value maps reveal structure on very small scales. Most of the 
surface density structure in these high resolution maps traces the dense gas structure in the 
cloud quite closely. However, as we increase n and approach measurements expected to be 
more robust in the face of Poisson statistics (for e xample, NN18 irnplies a surface density 
signal to noise ratio of 4 on any given measurement ( Casertano &: Hiitl 19851 )). most structure 
disappears, leaving behind only the two-peaked nature of the "double cluster" . 



•^In summary, maps are generated by measuring the area A of the smallest circle that contains the n nearest 
stars centered at regular grid intervals; the local surface density at the ith grid point is simply aj = (n— 1)/A i 
and the uncertainty in each measurement is equal to the density divided by \/n — 2 ( Casertano fc Hut IQSSh . 
We will be referencing nearest neigbor maps of differing n throughout this work, thus we adopt the notation 
NNn to represent this concept. 
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Fraction of Members with Disks 



It has been demonstrated previously that the fraction of cluster membe rs with disks 
roughly traces the mean ag e of the pre-main sequence cluster members (e.g. iHaisch et al. 



2OOII : [Hernandez et al.ll2007l ). Here we present a measurement of this quantity for NGC 1333 
in order to place it within the context of this trend. 

Based on the overall structure observed in the density maps, we have adopted a circular 
region labeled "Main" in Figure [H] that is clearly dominated by cluster members for mea- 
surement of the fraction of members that have disks. The circle is centered on the median 
Right Ascension and Declination of the Spitzer-identified YSOs, and has a radius of 5.5 ar- 
cminutes, or 0.4 pc. We have adopted a rectangular field star reference region, also labeled 
in Figure [TH where no YSOs are found and little structure is present in the extinction map 
or the SCUBA map. Finally, we have chosen several smaller regions of the cluster for similar 
analysis to examine if the disk fraction changes between the northern and southern peaks 
in stellar density, as well as the diffuse Class Il-dominated region to the east of the main 
cluster. 

Figure [15] demonstrates the Ks magnitude histogram (KMH) analysis that yields an 
estimate of the total membership, including diskless pre-main sequence stars, within the 
"Main" region down to the adopted magnitude limit of Ks = 14 for all sources that have 
at least H, Ks, [3.6], and [4.5] photor aetry. A detailed descrip tion of the application of this 
process can be found in Section 3 of iGutermuth et al.l (120051 ). In summary, we count the 
number of stars per 0.5 magnitude bin in the cluster region and in the reference region and 
take the difference, applying a correction to the reference region distribution to account for 
the mean difference in extinction and the ratio of the areas between the two regions. After 
field star contamination is accounted for, we estimate that 87 members are located within 
the "Main" cluster region. The a dopted Ks magnitu de limit, corresponding to a 1 Myr old, 
0.08 Mq pre-main sequence star (IBaraffe et al.lll998l ) behind Ak = 2 mag of extinction, was 
adopted to ensure uniform detection likelihood in sufficient near-IR and mid-IR wavelengths 
to enable confident classification of sources. Further, this cut, in combination with the 
i^-band completeness limit of 16.1, achieves an appr oximate mass- and extinction-limited 
sampling of cluster members ([Carpenter et al.l 119971 ). In this way, we effectively mitigate 
biases against detecting diskless members over members with disks (the longer wavelength 
Spitzer photometry is deeper that the near-IR photometry, enabling detection of both types 
of sources for a given near-IR magnitude cut) and detecting low- luminosity members with 
disks that may have been cut in the elimination of extragalactic contaminants (though the 
number of these sources appears to be small, see Fig. [6] in Section HTTj) . The total number 
of Spitzer-identified YSOs with Ks < 14 in the same region is 72, yielding a disk fraction of 
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83% ±11%. 

We have summarized the disk fraction measurements for all subregions in Table |H 
All measurements are consistent to within the uncertainty of the field star contamination 
estimates, suggesting that both components of the "double cluster" as wel l as the more 



diffuse population to the east are roughly coeval. For further comparison, IWilking et al. 



(I2OO4J ) estimated a disk fraction for the northern subcluster of 75% ± 20%, consistent with 
our results, using their K-band excess fraction with a correction factor applied to account 
for this method's incomplete census of stars with disks. 



7. ID Surface Density Profiles 

Most recent models of clustered star formation do not attempt to predict specific ra- 
dial density profile falloff behavior, as such models suggest that high dynamical activity 
dominates these environments, making any structure in a simulated cluster quite transient. 
However, recent N2H^ measurements of the radial velocitie s of cold cores in NGC 1333 show 
that the velocity dispersion is very low (jWalsh et al. 2006 ). An analysis of the distances of 



red MIPS sources in NGC 1333 and the rest of Perseus from submillimeter cores suggests 
that the protostellar so urces have a very low velocity dispersion relative to the dense gas 



( lJ0rgensen et al.ll2006bl ). Furthermore, the similarity in the nearest neighbor spacings be- 
tween the Class I and Class II sources is also indicative that the relative speeds of the cluster 
members is small and the Class II sources are particularly young. Therefore, the structure 
of NGC 1333 is particularly intriguing as it should offer a direct probe of an early phase of 
the clustered star formation process, before gravitational interactions have had a chance to 
erase structure in the member distribution. Here we measure the radial density profile for 
NGC 1333, as well as measure the degree of asymmetry in the cluster. 



7.1. Radial Density Profile: Annular Star Counts Method 



Unlike the regions examined in lGutermuth et al.l (120051 ). the density of cluster members 
relative to the field star density in the NGC 1333 region is too low to allow for any meaningful 
measurement of the radial profile of the cluster from the Ks-hand source list. Thus, we 
restrict this analysis of the character of the cluster radial profile to the one measured for the 
Spitzer-identified YSOs only. Given the large fraction of the total membership inferred to 
have disks, this should introduce minimal additional uncertainty to the analysis. 



Here we adopt the method of measuring the radial surface density profile described 
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in detail in iGutermuth et al.l (120051 ) . In summary, once a cluster center point is adopted, 
sources are counted in concentric annular regions of constant radial extent, and the counts 
are divided by the area of each annulus. The bin size used here is 1 arcminute, or roughly 
0.07 pc. However, choosing the center point is an uncertain task for clusters that are non- 
uniform and/or asymmetric. 

Given the degree of substructure in NGC 1333, we have measured radial profiles from 
four distinct center points: "Main" is the median of the Right Ascension and Declination 
of the YSOs; "Ring Center" is the visual center of the cluster, located at the center of the 
apparent ring traced by the dense gas and YSOs; "North" and "South" are the locations of 
peak density in the two main maxima of the double cluster. See Figure [16] for the resulting 
profiles. The YSOs in NGC 1333 are distributed roughly uniformly in the inner 0.3 pc 
radius region of the cluster. The choice of center point primarily affects the short radius 
behavior of the profile, depending on proximity to a local maximum. Also, the profiles made 
from offset center points ("North" and "South") have clear contamination from the other 
local maximum at ~0.55 pc. This ambiguity demonstrates that traditional radial profile 
techniques like this are not particularly well-suited for studies of young clusters. 



Gutermuth et al.l (120051 ) fit measured radial profiles with exponential decay functions, 
but noted without demonstration that the functional choice seemed to be arbitrary; one 
could fit any number of well-behaved functions to the cluster radial profiles measured for 
those clusters given the uncertainty introduced by field star contamination. The radial 
profiles derived from the NGC 1333 YSOs do not remotely resemble an exponential, nor 
even a single index power law, however. The dashed lines in Figure [16] are two power laws 
chosen to emulate the functional behavior of the "Main" and "Ring Center" profiles in the 
inner and outer radius regimes. Inside the knee at 0.3 pc, the profile is approximately fiat 
(power law index of a = 0). Outside that knee, the descent is precipitous, corresponding 
to a power law index of a = —3.3. Alternatively, one can also model the behavior of the 
profile with a slightly more complicated function, such as the red curve in Figure [16] This 
is a function of the form a = ctq^I + , with parameters ctq = 220 pc~^, tq = 0.3 pc, 

P = 6, and 7 = 3.3. 

The function above is very similar in functional form to the EFF p rofile often used t o 



characterize the radial profiles of young massive clusters in other galaxies (lElson et al.lll987l ). 
although the value of f3 in that work is kept constant at /? = 2. We allowed that value to 
vary, yielding a closer fit to the entire profile. For relatively extreme fitting parameters ( 
(To = 220 pc~^. To = 0.7 pc, arid 7 = 8.0), an EFF profile can be fit to these data. Similarly, 



we can fit a King profile (lKing]ll962l ) of central density k = 400 pc ^, core radius Vc = 0.4 pc. 



and truncation radius rt = 1.3 pc. See the right column of plots in Figure [16] for overlays of 
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these functions on the measured radial profiles. Regardless of the fitting results, the physical 
inferences one can normally derive from these models are not valid, as NGC 1333 is a very 
non-uniform, asymmetric system (see Section 17.21 below) that appears to be dominated by 
structure in the dense gas. 



7.2. Azimuthal Distribution Histogram 

We next apply another one dimensional analysis of structure, averaging rad ially rather 



than azimuthally in order to measure the azimuthal distribution histogram (ADH; iGutermuth et al 



20051 ) for the YSOs in NGC 1333. The ADH is presented in Figure [T71 Non-uniform az- 
imuthal structure is apparent, with peaks at PA~ 30° and ~ 210°. These peaks correspond 
to the "double cluster" local density maxima. To quantify the asymmetry we see here, we 
measure the azimuthal asymmetry parameter (AAP) of this ADH, defined as the ratio of the 



standard deviation of the histogram to the square root of its mean value in IGutermuth et al. 



( I2OO5I ). The AAP is meant to quantify the degree of asymmetry measured in an ADH with 
respect to the deviations expected from a normal, circularly symmetric cluster distribution, 
e.g. one characterized by Poisson variability only. Monte Carlo simulations presented in 



Gutermuth et al.l (120051 ) demonstrated that for clusters of ~ 100 sources, the behavior 
of the AAP of an unstructured cluster is reasonably approximated as a Gaussian with a 
mean value of ~0.9 and a standard deviation of ~0.2. Thus the AAP of 2 measured here 
for NGC 1333 is quite high, demonstrating that the cluster is significantly deviant from a 
circularly symmetric distribution. 



8. Discussion 

We have employed the mid-IR sensitivity of Spitzer to achieve a census of the YSO 
members of the NGC 1333 embedded cluster that was previously unattainable. Furthermore, 
we have shown that the sources identified in our Spitzer census represent a large fraction 
of the total cluster membership (83%). With this penetrating view of the cluster, we have 
performed several measurements of different aspects of the structure of the cluster that 
are not confused by field star contamination and are less biased by extinction effects than 
previous studies. 



We confirr a the double-peaked surfac e density morphology of the cluster reported in 
previous work (ILada. Alves. fc Ladalll996l ). with the caveat that this morphology is traced 
only by the more evolved Class II population. The two main density peaks are located in 
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local minima of the gas density distribution. In contrast, the protostars trace the dense gas 
in this region closely, as do a fraction of the Class II. That gas distribution appears to be 
a network of filamentary structures, connecting the two density peaks into a single, albeit 
complex, structure. Despite the difference in spatial distributions, however, the Class II and 
protostars have similarities in their nearest neighbor distance distributions, particularly in 
the location of the peak at spacings of 0.045 pc. These two evolutionary states are expected 
to differ in duration by an order of magnitude, thus a low overall velocity dispersion in 
the cluster and a very young Class II population are needed to account for the similarity 
remaining in the spacings among the two populations. 

A further implication of the lack of evidence for dynamical evolution of the cluster is the 
need for the dispersal of the molecular gas by the Class II sources on very short timescales. 
Given that the more evolved YSOs have not moved significantly from their birthsites, regions 
dominated by more evolved sources where gas densities are preferentially low, such as the 
double-peaks in the stellar distribution here, were once locations of high gas density. Thus 
it seems plausible that groups of low-mass stars are able to disrupt their dense natal gas 
locally. In this sense, we concur with previous work that has suggested that the activ e 



outflows in NGC 1333 are destroying the cloud (jSandell fc Kneell200ll : iQuillen et al.ll2005l ). 
With a virtual lack of low or high mass stars and dense gas in the center of the ring-like 
structure of the cloud and cluster core, we argue that this is in fact a primordial structure 
and not outflow-evacuated. 

Considering all sources regardless of evolutionary class, the cluster is clearly elongated, 
an expected result given the double-peaked nat ure of the stellar surface density distribution 



and embedded nature of the cluster as a whole (IGutermuth et al.ll2005l ). We have presented 
radial density profiles measured via several different methods, and all suggest a roughly 
uniform density distribution within a 0.3 pc radius, with a steep decline {a = —3.3) beyond. 
The fiat central distribution and sharp radial decline in the surface density profile at radii 
larger than 0.3 pc also suggests a rather limited amount of dynamical interaction in this 
cluster. If we analyze this within the framework of dynamics-generated structures like those 
of King or EFF models, we have to choose extreme fitting parameters just to get marginal 
agreement with the measured profiles. The radial profile we have measured here has two 
clear regimes were a simple power law matchs the density profile well, and the sharp knee 
transition between them is poorly matched by either the King or EFF profiles. Furthermore, 
given the asymmetry of the NGC 1333 cluster and the presense of a considerable amount 
of structured, dense gas, the dynamical information inferred from either of these models is 
unlikely to be accurate. Cloud geometry is the most tenable cause for the structure we have 
observed. As such, the relative motions of the stars must be fairly slow and the Class II 
population must be quite young in order to have preserved that structure into the Class II 
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evolutionary phase. 
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A. Extragalactic Contamination Characterization and Mitigation 



While IR-excess emission is a potent membership diagnostic for young and embed- 
ded clusters, there are many potential contaminants. For NGC 1333, and most clusters in 
the Spitzer Young Stel lar Clusters Surve y, the dominant source of contamination is from 
extragalactic source^. Stern et al. ( 2005 ) demonstrated that there were two classes of ex- 
tragalactic sources that have some form of IR-excess emission and were quite numerous once 
sufficiently deep integrations are obtained with IRAC. First, normal star-forming galaxies 
and narrow-line AGN have growing excess at 5.8 and 8.0 /xm, as PAH emission excited by 
young, high-mass stars dominates their largely stellar SEDs in these bandpasses. Second, 
broad-line AGN have red, non-stellar SEDs, resulting in IRAC colors that are very similar 
to those of bona fide YSOs. 



To aid in the development of mitigation strategies for the Spitzer Young Cluster survey 
and stu dies of nearby star formation in general, we have obtained the Bootes field IRAC data 
used in lStern et al.l (120051 ) from the Spitzer Science Center online public archive and reduced 
~7.7 square degrees of it using identical methods to those presented in Section [2l The total 



"^For regions closer to the galactic plane or nearer to the galactic center, evolved stars are also relatively 
bright and numerous, but this source of contamination is not expected to be significant for most nearby 
star-forming regions. Thus, wc do not addressed that in this paper. 
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exposure time in the mosaicked data is twice that of the NGC 1333 mosaics, and the field 
lacks significant structured nebulosity. Thus the Bootes data are marginally deeper than 
NGC 1333, moreso in the 5.8 and 8.0 fim mosaics than the 3.6 and 4.5 /xm ones. Without 
any attempt at cutting out contaminants, our IRAC four channel color-based classification 
scheme would sugges t that we had detect ed 12 68 Class I and 2641 Clas s II YSOs in this field; 
clearly a false result. iHarvey et al.l (120061 ) and I J0rgensen et al.l (l2006al ) suggested a basic cut 
developed for the c2d Legacy Survey to filter these sources: cut all sources with [4.5] < 14. 
This method is relatively effective for the Bootes data in that only 80 (10.4 deg~^) Class I 
and 108 (14 deg^^) Class II of the above numbers of IR-excess sources would bypass the 
filter. 

In exploring the efficacy of the c2d filtering method, we determined a few simple refine- 
ments that seem to improve the overall performance of the method. We present the modified 
method here, along with a characterization of the number density of sources per magnitude 
for verification of this and other extragalactic filtering methods. In Figures [T8ll20l we show 
select color-color and color-magnitude diagrams to give a sense of the regions of these dia- 
grams where these non-YSO IR-excess sources fall, their approximate number densities as 
a function of survey depth, and a suggested set of color and magnitude cuts that can filter 
out most of these sources while minimizing loss of bona fide YSOs. 



First, we exploit the result from IStern et al.l (120051 ) that galaxies dominated by PAH 
emission have mid-IR colors that occupy a relatively unique area of most IRAC color-color 
diagrams. One simply has to filter out those sources that fall within either of the following 
regions of two IRAC color-color spaces. For the [4.5] — [5.8] vs. [5.8] — [8.0] diagram (see the 
top plot in Figure [T8|) . we reject sources that follow all these conditions: 

[4.5]-[5.8]<^x([5.8]-[8.0]-l) 

[4.5] - [5.8] < 1.05 
[5.8] - [8.0] > 1 



For the [3.6] — [5.8] vs. [4.5] — [8.0] diagram (see bottom plot in Figure [T8|) . we reject 
sources that follow all these conditions: 

[3.6] - [5.8] < ^ X ([4.5] - [8.0] - 1) 

[3.6] - [5.8] < 1.5 
[4.5] - [8.0] > 1 
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Only a small number of true YSOs have been sh own to have significant PAH emission in 



other regions of young and forming groups of stars (IHartmann et al.ll2005l : iHernandez et al. 



20071 ). Since no PAH emission source in the Bootes field has [4.5] < 11.5, we require all 
sources removed to be dimmer than this value. This ensures that at least the PAH-rich 
YSOs known at this time will not be excluded using this cut. 



The broad-line AGN have mid-IR colors that are largely consistent with YSOs (jStern et al 



20051 ). Thus, to filter these sources from IRAC source lists for nearby star-forming regions, 
one must exploit the fact that the apparent magnitudes of the extragalactic sources are 
much dimmer than those of most YSOs within the nearest kiloparsec. By filtering out the 
PAH sources first, we can construct a filter that more closely traces the broad-line AGN 
distribution than the original c2d proposed cut. We have developed this modified cut after 
applying considerable effort to characterizing the contaminants fr om the Bootes field and 
finding ways to isolate them from YSOs. Following the c2d lead (IHarvey et al.l l2006l ) . we 
utilize the [4.5] vs. [4.5] — [8.0] color-magnitude diagram to flag likely AGN. However, with 
the PAH sources removed separately, we have refined the AGN cutoff, flagging all sources 
that follow all of the following three conditions: 



[4.5] - [8.0] > 0.5 

[4.5] > 13.5 + ([4.5] - [8.0] - 2.3)/0.4 
[4.5] > 13.5 

Additionally, a source flagged as a likely AGN must follow any one of the following three 
conditions: 



[4.5] > 14+ ([4.5] - [8.0] -0.5) 

[4.5] > 14.5 - ([4.5] - [8.0] - 1.2)/0.3 
[4.5] > 14.5 

In Figure [191 we demonstrate the results of the above cutoff for the Bootes data. This 
cut has been chosen to avoid cutting large sections of color or magnitude space where we 
expect YSOs in the nearest kiloparsec to fall, thus this method does allow some residual 
contamination. However, this set of cuts is a significant improvement over the c2d filtration 
method. When applied to the Bootes data, it flags 1485 sources as PAH emission sources and 
2766 as likely AGN, leaving behind 49 (6.4 deg'^) false Class O/I sources and 29 (3.8 deg"^) 
false Class II sources as residual contamination. 
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Since most of these contaminating sources (PAH and AGN) are relatively bright at 
8.0 yum, it appears that the band that limits sensitivity to these contaminant sources is the 
5.8 fim channel. Therefore, it is useful to characterize the [5.8] magnitude distribution of 
each class of source from the Bootes dataset to demonstrate the effectiveness of the method 
presented here and model the expected residual contamination in a statistical sense (see 
Figure [20]) . The differential surface density distributions in units of number of sources per 
square degree per half-magnitude bin are both represented reasonably well by the following 
power laws: 

apAH = 10-''-'^ X [5.8]i^-°2 

^ — in-20-90 y rp; 0119.63 

(^BL-AGN — -LU X [b.»J 

Also, while differential completeness limits are not very accurately determined using 
magnitude histograms such as these, the overall completeness "decay" behavior can be mea- 
sured from the PAH emission source distribution, since these sources are uniquely identifiable. 
This behavior can then be applied to the power-law model above for the more ambiguously 
identified AGN distribution in order to help to compute the total number of these sources 
down to the detection limit of a given set of photometry data. 
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Table 1. Summary of Photometric Data Parameters 



Bandpass 


Flux for Zero Mag.{Jy) 


Mag. for 90% Completeness'' 


Median Unc. at Completeness Limit 


Mag. for 1 DN/s^ 


Ax c 

Ak 


J 


1669 


16.8 


0.06 




2.50 


H 


980 


16.1 


0.08 




1.55 


Ks 


620 


15.3 


0.06 




1.00 


[3.6] 


280.9 


17.0 


0.08 


19.455 


0.63 


[4.5] 


179.7 


16.5 


0.09 


18.699 


0.53 


[5.8] 


115.0 


14.5 


0.11 


16.498 


0.49 


[8.0] 


64.13 


13.8 


0.12 


16.892 


0.49 


[24] 


7.3 


9.4 


0.08 


14.600 


0.43 



''2M ASS JHKs values inferred from 2MASS Online Documentation; Spitzer values computed using method detailed in ICutermuth et all 

liool). 

''Photometric calibration values account for adopted apertures. 
'^Adopted from iFlahertv et al. 



Table 2. Spitzer-identified YSOs: 2MASS, IRAC, and MIPS Photometry 







A a 


Decl.'' 


J 






PCs 


f3 fil 
[CS.DJ 




[O.OJ 


[8 nl 

[O.UJ 


[941 


1 


03 


28:39.09 


I O-L 


06:01.9 


... ± .. 




... + ... 


... + ... 


15.62 ± 0.04 


14.26 + 0.03 


13.31 + 0.05 


11.91 + 0.03 


6.17 + 0.01 


2 


03 


28:45.30 


I O-L 


05:42.1 


... ± .. 




16.75 + 0.15 


14.91 + 0.06 


13.17 + 0.02 


11.82 + 0.01 


11.37 + 0.02 


10.72 + 0.01 


3.73 + 0.01 


3 


03 


28:57.37 




14:16.2 


... ± .. 




... ± ... 


... ± ... 


... ± ... 


7.81 + 0.01 


6.60 + 0.01 


5.61 + 0.01 


0.47 + 0.01 


4 


03 


29:00.50 


4-31 


12:00.6 


... ± .. 




... ± ... 


... ± ... 


... ± ... 


15.11 + 0.09 


... ± ... 


... ± ... 


6.82 + 0.03 


5 


03 


29:04.06 


-ro± 


14:46.7 


... ± .. 




... ± ... 


... ± ... 


14.21 + 0.04 


12.81 + 0.04 


... ± ... 


11.54 + 0.03 


5.17 ± 0.06 


6 


03 


29:10.49 




13:30.8 


... ± .. 




... ± ... 


... ± ... 


... ± ... 


15.48 + 0.07 


... ± ... 


... ± ... 


6.13 ± 0.04 


7 


03 


29:11.25 


-1-31 


18:31.7 


... ± .. 




... ± ... 


... ± ... 


13.95 + 0.13 


11.61 + 0.06 


10.81 ± 0.06 


10.44 ± 0.06 


2.68 ± 0.02 


8 


03 


29:12.04 




13:02.0 


... ± .. 




... ± ... 


... ± ... 


13.95 ± 0.01 


10.47 + 0.01 


10.95 + 0.01 


10.70 + 0.04 


5.55 + 0.02 


9 


03 


29:13.60 


4-31 


13:58.2 


... ± .. 




... ± ... 


... ± ... 


15.69 + 0.04 


13.53 + 0.02 


12.83 + 0.06 


13.12 + 0.18 


6.00 + 0.03 


10 


03 


29:17.16 




27:46.1 


... ± .. 




... ± ... 


... ± ... 


14.28 ± 0.02 


12.78 + 0.01 


12.25 + 0.03 


12.16 + 0.05 


4.52 + 0.01 


11 


03 


29:23.48 


-f-O-L 


33:29.4 


... ± .. 




... ± ... 


... ± ... 


... ± ... 


12.35 + 0.01 


11.58 + 0.02 


10.90 ± 0.02 


3.75 ± 0.01 


12 


03 


28:32.55 


-1-31 


11:04.8 


... ± .. 




... ± ... 


14.49 + 0.12 


13.22 + 0.02 


12.31 + 0.01 


11.56 + 0.01 


10.64 ± 0.01 


5.29 + 0.01 


13 


03 


28:34.49 


-t-31 


00:51.2 


14.72 ± 


03 


13.09 ± 0.04 


11.89 + 0.03 


9.12 + 0.01 


8.20 + 0.01 


6.81 + 0.01 


6.20 + 0.01 


1.89 + 0.01 


14 


03 


28:34.52 


+ 31 


07:05.5 


... ± .. 




... ± ... 


... ± ... 


12.15 ± 0.01 


10.55 + 0.01 


9.37 + 0.01 


8.25 + 0.01 


5.18 + 0.01 


15 


03 


28:37.09 




13:31.1 


... ± .. 




... + ... 


12.05 ± 0.03 


9.81 + 0.01 


8.01 + 0.01 


6.67 + 0.01 


4.95 + 0.01 


... + ... 


16 


03 


28:39.69 


+ 31 


17:32.0 


... ± .. 




17.04 ± 0.20 


13.81 ± 0.03 


10.81 + 0.01 


9.69 ± 0.01 


8.88 + 0.01 


8.00 ± 0.01 


4.03 ± 0.01 


17 


03 


28:40.62 


+31 


17:56.5 


... ± .. 




... ± ... 


... ± ... 


16.93 + 0.18 


14.93 ± 0.05 


13.44 ± 0.05 


12.30 ± 0.05 


7.71 + 0.13 


18 


03 


28:43.27 


-roi 


17:33.1 


12.48 ± 


02 


10.84 ± 0.05 


9.51 + 0.04 


7.82 ± 0.01 


6.86 + 0.01 


5.91 + 0.01 


4.90 ± 0.01 


1.32 + 0.01 


19 


03 


28:48.76 


+^1 


16:08.9 


... ± .. 




16.84 + 0.16 


14.95 ± 0.05 


13.19 + 0.01 


12.39 + 0.01 


11.48 + 0.01 


10.42 + 0.01 


6.50 ± 0.02 


20 


03 


28:51.27 


+ 31 


17:39.5 


16.02 + 


05 


14.83 ± 0.08 


14.31 ± 0.06 


13.91 ± 0.03 


12.96 + 0.02 


12.09 + 0.02 


10.49 + 0.01 


6.02 ± 0.02 


21 


03 


28:55.52 


+ 31 


14:35.8 


... ± .. 




... ± ... 


... ± ... 


12.66 ± 0.01 


10.18 + 0.01 


8.96 + 0.01 


8.33 + 0.01 


0.35 + 0.01 


22 


03 


28:57.09 


+ '=11 


21:25.1 


... ± .. 




... ± ... 


16.37 + 0.14 


14.35 ± 0.03 


13.59 + 0.03 


12.87 + 0.04 


12.12 + 0.03 


... ± ... 


23 


03 


28:57.11 


+31 


19:11.9 


17.21 + 


08 


15.97 + 0.09 


15.24 + 0.05 


14.03 ± 0.04 


13.08 ± 0.04 


12.22 + 0.05 


11.12 + 0.05 


7.26 + 0.10 


24 


03 


28:57.99 


+31 


20:52.1 


... ± .. 




... ± ... 


... ± ... 


15.47 + 0.08 


14.37 ± 0.07 


13.09 ± 0.07 


11.89 + 0.04 


... ± ... 


25 


03 


28:58.42 


+31 


22:17.7 


17.99 ± 


16 


14.43 ± 0.04 


11.56 + 0.03 


9.60 + 0.01 


8.01 + 0.01 


7.04 ± 0.01 


5.96 ± 0.01 


2.28 + 0.01 


26 


03 


28:59.31 


+31 


20:08.2 


17.80 ± 


14 


... ± ... 


15.10 + 0.05 


14.21 ± 0.09 


13.25 + 0.07 


12.23 + 0.16 


11.18 + 0.09 


... ± ... 


27 


03 


29:01.54 


+31 


20:20.7 


... ± .. 




13.88 ± 0.08 


10.88 ± 0.04 


6.92 + 0.01 


5.73 + 0.01 


4.53 + 0.01 


3.50 + 0.01 


... ± ... 


28 


03 


29:03.32 


+31 


23:14.8 


17.63 ± 


11 


15.48 ± 0.06 


14.02 ± 0.04 


12.00 + 0.01 


10.98 + 0.01 


9.49 + 0.01 


7.93 + 0.01 


4.31 + 0.02 


29 


03 


29:03.76 


+31 


16:04.0 


11.74 + 


02 


9.74 ± 0.05 


8.15 + 0.02 


6.22 + 0.01 


5.02 + 0.01 


3.92 ± 0.01 


2.73 + 0.01 


... ± ... 


30 


03 


29:04.94 


+31 


20:38.6 


... ± .. 




14.22 ± 0.07 


12.65 + 0.04 


11.30 + 0.01 


10.52 + 0.01 


9.72 + 0.01 


8.78 + 0.01 


4.80 ± 0.04 


31 


03 


29:07.77 


+31 


21:57.3 


... ± .. 




13.80 ± 0.09 


10.43 ± 0.04 


6.65 + 0.01 


5.13 + 0.01 


4.13 + 0.01 


3.65 ± 0.03 


... ± ... 


32 


03 


29:08.95 


+31 


22:56.3 


16.45 ± 0.07 


13.76 ± 0.05 


11.78 ± 0.03 


9.98 ± 0.02 


9.05 ± 0.01 


8.63 + 0.11 


... ± ... 


... ± ... 


33 


03 


29:09.07 


+31 


21:29.1 


... ± .. 




16.10 + 0.10 


13.18 ± 0.03 


10.30 + 0.01 


8.39 + 0.01 


7.42 ± 0.01 


6.74 ± 0.03 


2.34 ± 0.05 


34 


03 


29:09.09 


+31 


23:05.7 


14.61 ± 


03 


12.90 ± 0.03 


11.84 + 0.03 


10.65 ± 0.02 


10.10 ± 0.02 


8.97 + 0.10 


7.33 + 0.14 


2.83 + 0.12 


35 


03 


29:10.71 


+31 


18:21.2 


... ± .. 




... ± ... 


... ± ... 


12.57 + 0.01 


10.82 + 0.01 


9.77 + 0.01 


9.18 + 0.01 


2.27 + 0.02 



Table 2 — Continued 



Index 


R.A.^ 


Decl.^ 


J 


H 


Ks 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 



36 


03; 


;29:12.95 


+31; 


;18:14.6 


... ± ... 




... + ... 


14.14 + 0.03 


9.28 + 0.01 


7.53 + 0.01 


6.64 + 0.01 


5.80 + 0.01 


2.60 + 0.01 


37 


03; 


;29:18.25 


+31; 


;23:19.9 


... ± ... 




... + ... 


... + ... 


13.53 + 0.05 


11.99 + 0.02 


10.63 + 0.02 


9.59 + 0.03 


4.17 + 0.02 


38 


03; 


;29:22.28 


+31; 


:13:54.5 


... ± ... 




... + ... 


... + ... 


14.26 ± 0.01 


13.29 + 0.01 


12.50 + 0.02 


11.14 + 0.01 


7.73 + 0.03 


39 


03; 


;29:31.54 


+31; 


:25:27.8 


17.18 ± 0.07 


16, 


.03 ± 0.08 


15.35 ± 0.06 


14.34 + 0.01 


13.56 + 0.01 


12.78 + 0.03 


12.17 + 0.06 


7.90 + 0.09 


40 


03; 


;28:15.19 


+31; 


:17:23.8 


16.51 ± 0.05 


15, 


.75 ± 0.07 


15.13 ± 0.05 


14.45 + 0.01 


13.98 + 0.01 


13.49 + 0.05 


12.82 + 0.05 


... + ... 


41 


03; 


;28:34.47 


+31; 


;17;43.2 


... ± ... 


17, 


.13 + 0.19 


15.20 + 0.05 


13.96 + 0.01 


13.39 + 0.01 


12.96 + 0.03 


12.14 + 0.03 


8.50 + 0.05 


42 


03; 


;28:34.85 


+31; 


;16:04.5 


17.55 + 0.10 


16, 


.37 + 0.10 


15.13 + 0.05 


13.69 ± 0.01 


13.10 + 0.01 


12.48 + 0.02 


11.61 + 0.02 


8.03 + 0.05 


43 


03; 


;28:38.76 


+31; 


;18:06.7 


... ± ... 


16, 


.83 + 0.15 


14.16 ± 0.03 


12.00 + 0.01 


10.97 + 0.01 


10.41 + 0.01 


9.56 + 0.01 


6.82 + 0.04 


44 


03; 


;28:43.23 


+31; 


;10:42.7 


15.29 + 0.03 


13 


.15 + 0.03 


12.04 + 0.02 


11.20 + 0.01 


10.88 + 0.01 


10.57 + 0.01 


10.23 + 0.01 


8.13 + 0.04 


45 


03; 


;28:43.56 


+31; 


17:36.5 


12.02 ± 0.03 


10, 


.79 ± 0.03 


9.85 + 0.05 


8.90 + 0.01 


8.38 + 0.01 


7.77 + 0.01 


6.64 + 0.01 


... + ... 


46 


03; 


;28:44.08 


+31; 


: 20:52. 9 


14.21 ± 0.02 


13, 


.21 ± 0.03 


12.61 ± 0.03 


11.90 + 0.01 


11.68 + 0.01 


11.27 + 0.01 


10.61 + 0.01 


8.07 + 0.03 


47 


03; 


;28:46.19 


+31; 


:16:38.7 


10.88 ± 0.02 


10, 


.01 ± 0.02 


9.69 + 0.02 


9.47 + 0.01 


9.17 + 0.01 


8.88 + 0.01 


8.30 + 0.01 


... + ... 


48 


03; 


;28:47.63 


+31; 


;24;06.1 


14.26 + 0.02 


12, 


.61 + 0.03 


11.64 + 0.02 


10.34 + 0.01 


9.97 + 0.01 


9.66 + 0.01 


9.27 + 0.01 


6.52 + 0.01 


49 


03; 


;28:47.82 


+31; 


;16:55.3 


12.93 + 0.02 


11 


.75 ± 0.03 


10.92 + 0.02 


9.94 + 0.01 


9.33 + 0.01 


9.01 + 0.01 


8.33 + 0.01 


5.60 + 0.02 


50 


03; 


;28:51.02 


+31; 


18:18.5 


11.36 ± 0.02 


10 


.07 + 0.02 


9.18 + 0.02 


8.-50 + 0.01 


7.84 + 0.01 


7.55 + 0.01 


6.49 + 0.01 


3.37 + 0.01 


51 


03; 


;28:51.07 


+31; 


16:32.6 


13.29 ± 0.02 


12, 


.50 ± 0.03 


12.14 + 0.02 


11.91 ± 0.01 


11.41 + 0.01 


11.07 + 0.01 


10.37 + 0.01 


8.29 + 0.07 


52 


03; 


;28:51.19 


+31; 


19:54.9 


11.72 ± 0.02 


10, 


.49 ± 0.02 


9.90 + 0.02 


9.08 + 0.01 


8.65 + 0.01 


8.36 + 0.01 


7.29 + 0.01 


4.31 + 0.01 


53 


03; 


;28:52.13 


+31; 


;15;47.2 


13.19 + 0.02 


12, 


.50 ± 0.03 


12.07 + 0.02 


11.48 + 0.01 


11.21+0.01 


10.77 + 0.01 


10.04 + 0.01 


7.79 + 0.07 


54 


03; 


;28:52.15 


+31; 


;22:45.4 


11.98 + 0.02 


11, 


.01 + 0.02 


10.56 ± 0.02 


10.18 + 0.01 


9.94 + 0.01 


9.56 + 0.01 


8.92 + 0.01 


6.32 + 0.02 


55 


03; 


;28:52.90 


+31; 


;16:26.6 


13.59 ± 0.02 


12, 


.87 + 0.03 


12.50 ± 0.02 


12.26 + 0.01 


11.85 + 0.01 


11.42 + 0.01 


10.82 + 0.02 


8.29 + 0.17 


56 


03; 


;28:53.58 


+31; 


;12:14.7 


... ± ... 


16 


.80 + 0.15 


15.38 + 0.06 


14.10 + 0.01 


13.58 + 0.01 


13.09 + 0.04 


12.34 + 0.04 


9.07 + 0.12 


57 


03; 


;28:53.93 


+31; 


:18:09.3 


14.57 + 0.03 


12, 


.22 ± 0.03 


10.92 ± 0.03 


10.09 + 0.01 


9.69 + 0.01 


9.15 + 0.01 


8.30 + 0.01 


4.82 + 0.01 


58 


03; 


;28:54.07 


+31; 


16:54.5 


13.02 ± 0.02 


12 


.05 ± 0.05 


11.61 ± 0.02 


11.19 ± 0.01 


10.65 + 0.01 


10.21 + 0.01 


9.31 + 0.01 


5.12 + 0.01 


59 


03; 


;28:54.61 


+31; 


16:51.3 


12.82 ± 0.02 


11 


.22 ± 0.03 


10.32 ± 0.02 


9.54 + 0.01 


8.79 + 0.01 


8.45 + 0.01 


7.61 + 0.01 


4.50 + 0.02 


60 


03; 


;28:54.92 


+31; 


;15;29.2 


16.02 ± 0.04 


14, 


.97 ± 0.05 


14.18 + 0.03 


13.55 + 0.01 


13.07 + 0.01 


12.65 + 0.04 


11.99 + 0.04 


... + ... 


61 


03; 


;28:55.07 


+31; 


;16:28.8 


13.46 ± 0.02 


11, 


.61 + 0.03 


10.51 ± 0.04 


9.23 + 0.01 


8.47 + 0.01 


8.14 + 0.01 


7.46 + 0.01 


3.91 + 0.02 


62 


03; 


;28:55.15 


+31; 


;16:24.8 


12.96 ± 0.03 


11, 


.59 ± 0.05 


10.81 ± 0.03 


9.93 + 0.01 


9.50 + 0.01 


9.08 + 0.01 


8.30 + 0.01 


... + ... 


63 


03; 


;28:56.09 


+31; 


;19:08.6 


... ± ... 




... + ... 


16.18 + 0.12 


11.22 + 0.01 


9.71 + 0.01 


8.84 + 0.01 


8.26 + 0.01 


6.14 + 0.04 


64 


03; 


;28:56.31 


+31; 


: 22:28.0 


15.53 ± 0.03 


13, 


.07 ± 0.03 


11.77 + 0.03 


10.61 ± 0.01 


10.00 + 0.01 


9.27 + 0.01 


8.36 + 0.01 


4.75 + 0.03 


65 


03; 


;28:56.64 


+31; 


18:35.7 


12.30 + 0.02 


10 


.69 ± 0.02 


9.70 + 0.02 


8.83 + 0.01 


8.35 + 0.01 


8.25 + 0.01 


7.55 + 0.01 


3.67 + 0.01 


66 


03; 


;28:56.93 


+31; 


: 20:48. 6 


15.47 + 0.03 


14, 


.48 ± 0.04 


13.87 + 0.03 


13.31 ± 0.01 


13.17 + 0.01 


12.77 + 0.03 


11.89 + 0.04 


... + ... 


67 


03; 


;28:56.95 


+31; 


;16;22.3 


13.72 + 0.02 


11, 


.89 ± 0.03 


11.03 + 0.02 


10.61 + 0.01 


10.31+0.01 


9.92 + 0.01 


8.93 + 0.01 


3.97 + 0.02 


68 


03; 


;28:57.17 


+31; 


;15:34.6 


15.38 + 0.03 


14 


.00 ± 0.03 


13.20 ± 0.03 


12.61 + 0.01 


11.90 + 0.01 


11.26 + 0.02 


10.27 + 0.01 


7.17 + 0.09 


69 


03; 


;28:57.70 


+31; 


19:48.1 


13.05 ± 0.02 


11 


.95 ± 0.03 


11. ,38 ± 0.03 


10.62 ± 0.01 


10.10 + 0.01 


9.71 + 0.01 


8.67 + 0.01 


5.76 + 0.05 


70 


03; 


;28:58.24 


+31; 


:22:09.4 


16.04 ± 0.04 


14, 


.32 ± 0.04 


13.32 ± 0.03 


12.46 ± 0.02 


11.96 + 0.03 


11.36 + 0.07 


... + ... 


... + ... 



Table 2 — Continued 



Index 




R.A.* 


Decl.=- 


J 


H 


Ks 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24] 




71 


03; 


;28; 


;58.25 


+31; 


:22; 


:02.1 


14.55 ± 0.02 


13.09 ± 0.03 


12.22 + 0.02 


11.30 + 0.01 


10.84 + 0.01 


10.30 + 0.01 


9.59 + 0.04 


... + .. 




72 


03; 


;28; 


;58.43 


+31; 


:22: 


:56.8 


15.31 + 0.03 


14.24 ± 0.04 


13.65 + 0.03 


13.00 + 0.01 


12.73 + 0.01 


12.58 + 0.05 


12.07 + 0.04 


... + .. 




73 


03; 


;28; 


;59.32 


+31: 


;15: 


:48.5 


16.16 ± 0.04 


12.54 + 0.03 


10.56 ± 0.02 


8.69 + 0.01 


8.01 + 0.01 


7.50 + 0.01 


6.40 ± 0.01 


2.46 + 0.01 


74 


03; 


;28; 


;59.56 


+31: 


;21: 


:46.8 


12.37 + 0.02 


10.96 ± 0.03 


9.98 + 0.05 


8.85 + 0.01 


8.40 + 0.01 


8.14 + 0.01 


7.51 + 0.01 


4.77 + 0.03 


75 


03; 


;29; 


;00.15 


+31: 


;21: 


;09.3 


16.61 ± 0.05 


14.60 + 0.04 


13.35 ± 0.03 


12.11 ± 0.01 


11.68 + 0.01 


11.41 + 0.02 


10.62 + 0.03 


6.46 + 0.12 


76 


03; 


;29; 


;02.16 


+31: 


;16: 


;11.4 


14.51 + 0.02 


13.87 + 0.03 


13.62 ± 0.03 


13.33 ± 0.04 


12.76 + 0.06 


... + ... 


... + ... 


... + .. 




77 


03; 


;29; 


;02.79 


+31; 


:22; 


:17.4 


17.07 + 0.07 


14.88 ± 0.04 


13.54 ± 0.03 


12.20 + 0.01 


11.56 + 0.01 


11.07 + 0.02 


10.46 + 0.03 


... + .. 




78 


03; 


;29; 


;03.13 


+31; 


:22; 


:38.2 


14.05 ± 0.02 


12.78 ± 0.03 


11.73 + 0.02 


10.25 + 0.01 


9.60 + 0.01 


9.21 + 0.01 


8.23 + 0.01 


5.00 + 0.07 


79 


03; 


;29; 


;03.20 


+31: 


:25: 


:45.3 


15.66 + 0.03 


14.62 ± 0.04 


13.98 ± 0.03 


13.22 + 0.01 


12.64 + 0.01 


12.24 + 0.03 


11.56 + 0.05 


7.34 + 0, 


,09 


80 


03; 


;29; 


;03.39 


+31: 


;18: 


;40.1 


15.86 ± 0.03 


14.65 ± 0.04 


13.97 + 0.03 


13.26 ± 0.01 


13.02 + 0.02 


12.68 + 0.04 


12.37 + 0.04 


... + .. 




81 


03; 


;29; 


;03.42 


+31: 


;25: 


;14.5 


16.32 ± 0.05 


14.87 ± 0.05 


13.93 ± 0.04 


13.04 + 0.01 


12.67 + 0.01 


12.15 + 0.02 


11.50 + 0.04 


... + .. 




82 


03; 


;29; 


;03.86 


+31: 


:21: 


;48.8 


11.43 ± 0.02 


10.16 ± 0.03 


9.36 + 0.04 


8.05 + 0.01 


7.62 + 0.01 


7.24 + 0.01 


6.32 + 0.01 


3.50 + 0, 


,03 


83 


03; 


;29; 


;03.93 


+31; 


:23: 


:31.2 


17.20 + 0.08 


16.06 + 0.07 


15.17 + 0.06 


13.95 + 0.01 


13.45 + 0.01 


12.90 + 0.07 


12.42 + 0.17 


... + .. 




84 


03; 


;29; 


; 04.67 


+31: 


:16: 


:59.2 


15.61 + 0.03 


13.95 + 0.03 


12.75 + 0.02 


11.40 + 0.01 


10.85 + 0.01 


10.24 + 0.01 


9.22 + 0.01 


5.64 + 0, 


,14 


85 


03; 


;29; 


;04.73 


+31: 


;11: 


;35.0 


... ± ... 


... + ... 


14.31 ± 0.03 


12.04 + 0.01 


11.03 + 0.01 


10.41 + 0.01 


9.82 + 0.01 


7.39 + 0, 


,04 


86 


03; 


;29; 


;05.64 


+31: 


;20: 


;10.7 


17.20 ± 0.07 


16.34 + 0.09 


15.65 ± 0.08 


13.98 ± 0.02 


13.39 + 0.02 


12.48 + 0.02 


11.75 + 0.03 


... + .. 




87 


03; 


;29; 


;05.67 


+31: 


;21: 


;33.8 


... ± ... 


15.25 ± 0.05 


13.31 ± 0.03 


11.70 + 0.01 


10.92 + 0.01 


10.33 + 0.02 


9.75 + 0.05 


... + .. 




88 


03; 


;29; 


;05.76 


+31: 


;16: 


;39.7 


14.29 + 0.02 


11.58 + 0.03 


9.82 + 0.06 


8.56 + 0.01 


7.82 + 0.01 


7.16 + 0.01 


6.32 + 0.01 


3.71 + 0, 


,05 


89 


03; 


;29:06.32 


+31; 


:13:46.5 


... ± ... 


14.93 ± 0.04 


12.75 + 0.02 


11.06 + 0.01 


10.29 + 0.01 


9.70 + 0.01 


9.00 + 0.01 


5.97 + 0.03 


90 


03; 


;29; 


;06.92 


+31; 


:29; 


:57.1 


16.70 ± 0.05 


15.55 + 0.05 


14.83 ± 0.04 


14.06 ± 0.01 


13.79 + 0.01 


13.48 + 0.04 


12.93 + 0.05 


9.76 + 0, 


,11 


91 


03; 


;29; 


;07.94 


+31: 


:22: 


:51.6 


12.98 + 0.02 


11.17 + 0.03 


10.13 + 0.04 


8.93 + 0.01 


8.39 + 0.01 


7.95 + 0.01 


7.38 + 0.05 


... + .. 




92 


03; 


;29; 


;08.95 


+31: 


;26: 


:24.1 


17.14 + 0.07 


15.99 ± 0.07 


15.73 ± 0.08 


13.38 ± 0.01 


12.65 + 0.01 


12.07 + 0.02 


11.26 + 0.03 


7.34 + 0, 


,10 


93 


03; 


;29; 


;09.33 


+31: 


;21: 


;04.2 


16.38 ± 0.04 


14.34 ± 0.03 


13.19 ± 0.03 


12.28 ± 0.01 


11.92 + 0.01 


11.36 + 0.07 


10.20 + 0.07 


... + .. 




94 


03; 


;29; 


;09.47 


+31: 


;27: 


;21.0 


14.11 ± 0.02 


13.18 ± 0.03 


12.67 + 0.03 


12.06 ± 0.01 


11.71 + 0.01 


11.34 + 0.01 


10.41 + 0.01 


7.82 + 0, 


,12 


95 


03; 


;29; 


;09.63 


+31: 


;22: 


;56.5 


11.26 + 0.02 


10.14 + 0.05 


9.50 + 0.03 


9.13 + 0.01 


8.92 + 0.01 


8.33 + 0.08 


7.19 + 0.19 


... + .. 




96 


03; 


;29; 


;10.18 


+31; 


:27; 


:15.9 


15.51 + 0.03 


14.69 ± 0.04 


14.20 ± 0.03 


13.74 ± 0.01 


13.54 + 0.01 


13.50 + 0.06 


... + ... 


... + .. 




97 


03; 


;29; 


;10.40 


+31; 


:21; 


:59.3 


9.34 + 0.02 


7.95 + 0.04 


7.17 + 0.02 


6.70 + 0.01 


6.54 + 0.01 


5.95 + 0.09 


4.63 + 0.14 


... + .. 




98 


03; 


;29; 


;10.46 


+31: 


:23: 


:35.0 


15.65 + 0.03 


13.85 ± 0.03 


12.85 + 0.03 


11.87 + 0.01 


11.41+0.01 


10.90 + 0.03 


10.30 + 0.10 


... + .. 




99 


03; 


;29; 


;10.82 


+31: 


:16: 


:42.7 


15.97 + 0.03 


14.47 + 0.04 


13.20 ± 0.02 


11.60 ± 0.01 


10.90 + 0.01 


10.16 + 0.01 


9.05 + 0.01 


5.87 + 0, 


,03 


100 


03; 


;29; 


;11.33 


+31: 


;22: 


;57.1 


14.84 + 0.02 


13.46 ± 0.03 


12.58 ± 0.02 


11.50 + 0.01 


11.04 + 0.01 


10.72 + 0.12 


... + ... 


... + .. 




101 


03; 


;29; 


;11.64 


+31: 


;20: 


;37.7 


15.34 + 0.03 


13.60 ± 0.03 


12.69 ± 0.03 


11.93 ± 0.01 


11.60 ± 0.01 


11.16 + 0.02 


10.53 + 0.03 


6.02 + 0, 


,08 


102 


03; 


;29; 


;11.77 


+31: 


;26: 


;09.7 


16.87 + 0.06 


15.77 + 0.06 


14.76 ± 0.04 


13.85 + 0.01 


13.29 + 0.01 


12.64 + 0.07 


11.22 + 0.08 


6.35 + 0, 


,12 


103 


03; 


;29; 


;11.87 


+31: 


:21: 


:27.1 


... ± ... 


15.52 + 0.05 


12.69 + 0.03 


10.58 + 0.01 


9.69 + 0.01 


9.00 + 0.03 


8.25 + 0.06 


... + .. 




104 


03; 


;29; 


;12.90 


+31: 


:23: 


;29.5 


13.43 ± 0.02 


12.58 ± 0.03 


12.10 ± 0.03 


11.29 ± 0.01 


10.78 + 0.01 


10.62 + 0.03 


10.09 + 0.07 


... + .. 




105 


03; 


;29; 


;13.04 


+31: 


;17: 


;38.4 


15.22 ± 0.02 


14.59 ± 0.04 


14.19 ± 0.03 


13.55 ± 0.01 


13.10 + 0.01 


12.76 + 0.05 


12.08 + 0.04 


... + .. 





Table 2 — Continued 



Index 




Decl.=- 


J 


n 




[3.bJ 


[4.5J 


[5.SJ 


[8.UJ 






106 


03:29:13.13 


+31: 


;22:52.9 


13.04 ± 0.02 


11.24 + 0.03 


10.09 ± 0.05 


8.83 ± 0.01 


8.38 ± 0.01 


7.95 + 0.01 


7.50 + 0.01 


5.61 + 0, 


,19 


107 


03:29:15.34 


+31: 


;29:34.6 


... ± ... 


... + ... 


15.83 + 0.09 


14.00 + 0.01 


13.39 + 0.01 


12.66 + 0.03 


11.90 + 0.03 


9.21 + 0. 


,08 


108 


03:29:15.66 


+31: 


:19:11.1 


... ± ... 


... + ... 


16.47 + 0.15 


14.99 + 0.02 


14.56 ± 0.02 


14.08 ± 0.08 


13.47 ± 0.09 


... + .. 




109 


03:29:16.59 


+31: 


: 23:49. 6 


13.17 ± 0.02 


11.82 ± 0.03 


11.16 ± 0.02 


10.32 ± 0.01 


9.73 ± 0.01 


9.21 + 0.01 


8.45 + 0.04 


... + .. 




110 


03:29:16.81 


+31: 


:23:25.4 


15.25 ± 0.03 


14.10 ± 0.03 


13.52 ± 0.03 


12.85 ± 0.01 


12.42 ± 0.01 


11.79 ± 0.02 


10.82 ± 0.02 


6.36 + 0, 


,10 


111 


03:29:17.66 


+31: 


;22:45.2 


9.95 ± 0.02 


8.94 + 0.02 


8.35 ± 0.02 


7.29 + 0.01 


6.82 ± 0.01 


6.28 + 0.01 


5.54 + 0.01 


2.89 + 0, 


,01 


112 


03:29:17.76 


+31: 


;19:48.2 


14.78 ± 0.03 


13.65 ± 0.03 


12.97 + 0.03 


12.24 + 0.01 


11.87 + 0.01 


11.53 + 0.01 


10.93 ± 0.01 


8.09 + 0, 


,05 


113 


03:29:18.65 


+31: 


;20:17.9 


... ± ... 


... ± ... 


14.61 ± 0.04 


13.55 + 0.01 


12.89 + 0.01 


12.27 + 0.02 


11.06 + 0.01 


6.30 + 0. 


,03 


114 


03:29:18.72 


+31: 


;23:25.5 


11.37 + 0.02 


10.38 ± 0.09 


10.26 ± 0.02 


9.58 + 0.01 


9.10 + 0.01 


8.35 + 0.01 


7.06 + 0.01 


3.81 + 0, 


,04 


115 


03:29:20.05 


+31: 


:24:07.6 


... ± ... 


14.69 ± 0.05 


12.16 ± 0.03 


9.97 ± 0.01 


8.94 + 0.01 


8.32 + 0.04 


7.32 + 0.08 


3.22 + 0, 


,07 


116 


03:29:20.42 


+31: 


:18:34.3 


14.03 ± 0.02 


11.79 ± 0.03 


10.39 ± 0.02 


8.57 ± 0.01 


7.95 ± 0.01 


7.49 + 0.01 


6.59 + 0.01 


3.08 + 0, 


,01 


117 


03:29:21.30 


+31: 


;23:46.5 


17.48 ± 0.10 


14.66 ± 0.04 


13.07 ± 0.03 


11.68 ± 0.01 


11.19 ± 0.01 


10.72 ± 0.04 


10.09 ± 0.12 


... + .. 




118 


03:29:21.55 


+31: 


:21:10.5 


12.41 + 0.02 


11.72 + 0.03 


11.43 + 0.02 


11.01 + 0.01 


10.68 + 0.01 


10.32 + 0.01 


9.62 + 0.01 


6.74 + 0, 


,03 


119 


03:29:21.87 


+31: 


;15:36.4 


11.10 + 0.02 


10.05 + 0.03 


9.45 ± 0.03 


8.59 + 0.01 


8.20 + 0.01 


7.89 + 0.01 


7.41 +0.01 


4.82 + 0, 


,01 


120 


03:29:23.15 


+31: 


: 20:30. 5 


12.35 ± 0.02 


11.64 + 0.03 


11.30 ± 0.02 


10.51 ± 0.01 


10.16 ± 0.01 


9.77 + 0.01 


8.93 + 0.01 


5.92 + 0, 


,01 


121 


03:29:23.23 


+31: 


: 26:53.1 


13.52 ± 0.02 


12.67 ± 0.03 


12.24 ± 0.03 


11.63 ± 0.01 


11.09 ± 0.01 


10.72 ± 0.01 


10.02 ± 0.01 


7.51 + 0, 


,14 


122 


03:29:24.08 


+31: 


:19:57.8 


15.55 + 0.04 


14.25 + 0.04 


13.59 + 0.03 


12.87 + 0.01 


12.44 + 0.01 


11.88 ± 0.01 


10.75 ± 0.01 


5.03 + 0, 


,01 


123 


03:29:25.92 


+31: 


;26:40.1 


11.13 ± 0.03 


10.11 ± 0.04 


... ± ... 


8.92 + 0.01 


8.14 + 0.01 


7.67 + 0.01 


6.32 + 0.01 


4.31 + 0, 


,01 


124 


03:29:28.01 


+31: 


;25:11.0 


17.42 ± 0.09 


15.94 ± 0.07 


14.89 ± 0.04 


13.87 + 0.01 


13.37 + 0.01 


12.85 + 0.05 


12.28 + 0.11 


... + .. 




125 


03:29:29.79 


+31: 


;21:02.8 


12.66 ± 0.02 


11.63 + 0.03 


11.18 + 0.02 


10.69 ± 0.01 


10.43 ± 0.01 


9.92 + 0.01 


8.95 + 0.01 


6.08 + 0.03 


126 


03:29:30.39 


+31: 


;19:03.4 


12.12 + 0.02 


11.37 + 0.03 


11.06 + 0.02 


10.49 ± 0.01 


10.12 + 0.01 


9.69 + 0.01 


8.93 + 0.01 


6.04 + 0, 


,01 


127 


03:29:32.56 


+31: 


: 24:37.0 


14.31 + 0.02 


12.90 + 0.03 


11.85 + 0.02 


10.52 ± 0.01 


10.06 ± 0.01 


9.80 + 0.01 


9.26 + 0.01 


6.54 + 0, 


,05 


128 


03:29:32.86 


+31: 


: 27:12. 7 


13.37 ± 0.02 


12.66 + 0.03 


12.32 + 0.03 


11.76 ± 0.01 


11.46 ± 0.01 


11.23 ± 0.01 


10.73 ± 0.02 


8.51 + 0, 


,10 


129 


03:29:35.70 


+31: 


:21:08.6 


... ± ... 


... ± ... 


... ± ... 


15.06 ± 0.02 


13.96 ± 0.01 


13.33 ± 0.05 


12.77 + 0.05 


... + .. 




130 


03:29:37.63 


+31: 


;02:49.3 


14.68 ± 0.02 


14.13 + 0.03 


13.77 + 0.03 


13.18 + 0.01 


13.01 + 0.01 


12.48 + 0.03 


11.76 + 0.03 


9.03 + 0, 


,05 


131 


03:29:37.72 


+31: 


;22:02.6 


13.95 ± 0.02 


13.37 ± 0.02 


12.97 + 0.02 


12.35 + 0.01 


12.06 + 0.01 


11.71 + 0.01 


11.00 + 0.01 


7.34 + 0, 


,03 


132 


03:29:44.15 


+31: 


; 19:47.5 


17.46 + 0.10 


16.59 + 0.12 


15.16 + 0.05 


13.38 ± 0.01 


12.67 + 0.01 


12.07 + 0.01 


11.46 ± 0.02 


9.02 + 0, 


,06 


133 


03:29:54.03 


+31: 


:20:53.1 


11.93 + 0.02 


11.13 + 0.02 


10.60 ± 0.02 


9.90 + 0.01 


9.50 + 0.01 


9.15 + 0.01 


8.27 + 0.01 


4.46 + 0, 


,01 


134 


03:29:02.69 


+31: 


:19:05.8 


17.71 + 0.12 


15.31 + 0.05 


13.63 ± 0.03 


12.21 ± 0.01 


11.50 ± 0.01 


11.30 ± 0.01 


11.17 + 0.02 


7.39 + 0, 


,06 


135 


03:29:09.41 


+31: 


: 14: 14.1 


... ± ... 


... ± ... 


... ± ... 


12.44 + 0.01 


10.91 ± 0.01 


10.26 ± 0.01 


10.16 ± 0.04 


7.43 + 0, 


,10 


136 


03:29:26.79 


+31: 


: 26:47.7 


10.83 + 0.02 


9.99 ± 0.03 


9.68 ± 0.03 


9.58 ± 0.01 


9.40 ± 0.01 


9.28 + 0.01 


8.96 + 0.01 


5.91 + 0, 


,02 


137 


03:29:29.26 


+31: 


;18:34.8 


12.58 + 0.02 


11.40 + 0.03 


10.98 ± 0.02 


10.67 + 0.01 


10.65 + 0.01 


10.52 + 0.01 


10.05 + 0.01 


4.71 + 0, 


,01 



''J2000 coordinates 



Table 3. Spitzer-identi&ed YSOs: Addendum 
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R.A/ 


Decl.^ 


Lit. Matches 


Class'" Ak"" 
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03:28:58.42 


+31:22:17.7 


LAL 166 




1.24 


26 


03:28:59.31 


+31:20:08.2 


ASR 116; LAL 170 




0.65 


27 


03:29:01.54 


+31:20:20.7 


SVS 12; ASR 114 




1.11 


28 


03:29:03.32 


+31:23:14.8 


LAL 191 




1.94 


29 


03:29:03.76 


+31:16:04.0 


SVS 13; ASR 1; LAL 196 


I 1.01 


1.16 


30 


03:29:04.94 


+31:20:38.6 


LAL 204 




0.06 


31 


03:29:07.77 


+31:21:57.3 


LAL 213 




0.56 


32 


03:29:08.95 


+31:22:56.3 


LAL 221 


I 1.79 




33 


03:29:09.07 


+31:21:29.1 


LAL 223 




1.12 


34 


03:29:09.09 


+31:23:05.7 


LAL 222 


I 0.91 


1.07 


35 


03:29:10.71 


+31:18:21.2 






0.97 



Table 3 — Continued 



Index R.A.^ Decl.^ Lit. Matches Class"^ Ak" aiRAc"^ 



36 


03 


29 


12 


95 


+31 


18 


14 


6 


ASR 30; LAL 261 


I 




1.03 


37 


03 


29 


18 


25 


+31 
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1.69 
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0.69 
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31 
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15 


19 
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17 


23 
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II 


0.04 


-0.96 


41 


03 


28 


34 


47 


+31 


17 


43 


2 
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-0.78 


42 


03 


28 


34 


85 


+31 


16 


04 


5 




II 




-0.44 


43 


03 


28 


38 


76 


+31 


18 


06 


7 


LAL 63 


II 




-0.12 


44 


03 


28 


43 


23 


+31 


10 


42 


7 


LAL 80 


II 


1.51 


-1.72 


45 


03 


28 


43 


56 


+31 


17 


36 


5 


ASR 127 


II 


0.41 


-0.24 


46 


03 


28 


44 


08 


+31 


20 


52 


9 


LAL 83 


II 


0.22 


-1.32 


47 


03 


28 


46 


19 


+31 


16 


38 


7 


ASR 128; LAL 89 


II 


0.27 


-1.51 


48 


03 


28 


47 


63 


+31 


24 


06 


1 




II 


0.83 


-1.61 


49 


03 


28 


47 


82 


+31 


16 


55 


3 


SVS 17; ASR 111; LAL 97 


II 


0.27 


-1.05 


50 


03 


28 


51 


02 


+31 


18 


18 


5 


SVS 10; ASR 122; LAL 106 


II 


0.41 


-0.62 


51 


03 


28 


51 


07 


+31 


16 


32 


6 


ASR 44; LAL 107 


II 


0.04 


-1.09 


52 


03 


28 


51 


19 


+31 


19 


54 


9 


ASR 125; LAL 110 


II 


0.57 


-0.83 


53 


03 


28 


52 


13 


+31 


15 


47 


2 


ASR 45; LAL 125 


II 


0.00 


-1.17 


54 


03 


28 


52 


15 


+31 


22 


45 


4 


LAL 120 


II 


0.30 


-1.36 


55 


03 


28 


52 


90 


+31 


16 


26 


6 


ASR 46; LAL 128 


II 


0.00 


-1.18 


56 


03 


28 


53 


58 


+31 


12 


14 


7 


ASR 97 


II 




-0.82 


57 


03 


28 


53 


93 


+31 


18 


09 


3 


ASR 40; LAL 129 


II 


1.95 


-0.76 


58 


03 


28 


54 


07 


+31 


16 


54 


5 


ASR 42; LAL 131 


II 


0.33 


-0.70 


59 


03 


28 


54 


61 


+31 


16 


51 


3 


SVS 18; ASR 43; LAL 136 


II 


0.93 


-0.70 


60 


03 


28 


54 


92 


+31 


15 


29 


2 


ASR 109; LAL 138 


II 


0.19 


-1.06 


61 


03 


28 


55 


07 


+31 


16 


28 


8 


ASR 107; LAL 141 


II 


1.20 


-0.89 


62 


03 


28 


55 
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+31 


16 


24 


8 


ASR 108 


II 


0.65 


-0.97 


63 


03 


28 


56 


09 
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19 


08 
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II 




0.47 


64 


03 


28 


56 


31 


+31 


22 


28 





LAL 147 


II 


2.19 


-0.22 


65 


03 


28 


56 


64 


+31 


18 


35 


7 


ASR 120; LAL 150 


II 


0.80 


-1.45 


66 


03 


28 


56 


93 


+31 


20 


48 


6 


LAL 152 


II 


0.05 


-1.17 


67 


03 


28 


56 


95 


+31 


16 


22 


3 


SVS 15; ASR 118; LAL 154 


II 


1.30 


-0.91 


68 


03 


28 


57 


17 


+31 


15 


34 


6 


ASR 17; LAL 156 


II 


0.70 


-0.16 


69 


03 


28 


57 


70 


+31 


19 


48 


1 


ASR 113 


II 


0.32 


-0.63 


70 


03 


28 


58 


24 


+31 


22 


09 


4 


LAL 164 


II 


0.67 





Table 3 — Continued 



Index 




ueci. 


Lit. Matches 




OlIRAC 


71 


03:28:58.25 


+31:22:02.1 


LAL 163 


II 0.43 


-0.85 


72 


03:28:58.43 


+31:22:56.8 


LAL 165 


II 0.69 


-1.80 


73 


03:28:59.32 


+31:15:48.5 


SVS 16; ASR 106; LAL 171 


II 3.53 


-0.25 


74 


03:28:59.56 


+31:21:46.8 


LAL 173 


II 0.44 


-1.34 


75 


03:29:00.15 


+31:21:09.3 


LAL 175 


II 1.01 


-1.16 


76 


03:29:02.16 


+31:16:11.4 


ASR 3; LAL 184 


II 0.00 




77 


03:29:02.79 


+31:22:17.4 


LAL 186 


II 1.20 


-0.87 


78 


03:29:03.13 


+31:22:38.2 


LAL 189 


II 0.37 


-0.57 


79 


03:29:03.20 


+31:25:45.3 


LAL 190 


II 0.31 


-0.96 


80 


03:29:03.39 


+31:18:40.1 


ASR 63; LAL 193 


II 0.71 


-1.79 


81 


03:29:03.42 


+31:25:14.5 


LAL 192 


II 0.94 


-1.04 


82 


03:29:03.86 


+31:21:48.8 


LAL 195 


II 0.54 


-0.87 


83 


03:29:03.93 


+31:23:31.2 






-1.06 


84 


03:29:04.67 


+31:16:59.2 


ASR 105; LAL 203 


II 0.63 


-0.32 


85 


03:29:04.73 


+31:11:35.0 


ASR 99; LAL 205 




-0.34 


86 


03:29:05.64 


+31:20:10.7 




II 


-0.20 


87 


03:29:05.67 


+31:21:33.8 


LAL 206 


II 


-0.61 


88 


03:29:05.76 


+31:16:39.7 


SVS 14; ASR 7; LAL 207 


II 2.19 


-0.28 


89 


03:29:06.32 


+31:13:46.5 


ASR 53; LAL 208 




-0.50 


90 


03:29:06.92 


+31:29:57.1 




II 0.35 


-1.54 


91 


03:29:07.94 


+31:22:51.6 


LAL 215 


II 1.16 


-1.07 


92 


03:29:08.95 


+31:26:24.1 


LAL 220 


II 0.57 


-0.42 


93 


03:29:09.33 


+31:21:04.2 


LAL 225 


II 1.51 


-0.43 


94 


03:29:09.47 


+31:27:21.0 




II 0.26 


-0.96 


95 


03:29:09.63 


+31:22:56.5 


SVS 7; LAL 226 


II 0.35 


-0.56 


96 


03:29:10.18 


+31:27:15.9 


LAL 228 


II 0.00 




97 


03:29:10.40 


+31:21:59.3 


SVS 3; LAL 230 


II 0.62 


-0.40 


98 


03:29:10.46 


+31:23:35.0 


LAL 231 


II 1.07 


-1.02 


99 


03:29:10.82 


+31:16:42.7 


ASR 23; LAL 235 


II 0.64 


0.09 


100 


03:29:11.33 


+31:22:57.1 


LAL 239 


II 0.66 




101 


03:29:11.64 


+31:20:37.7 


ASR 85; LAL 245 


II 1.27 


-1.20 


102 


03:29:11.77 


+31:26:09.7 


LAL 246 




0.16 


103 


03:29:11.87 


+31:21:27.1 


LAL 248 




-0.18 


104 


03:29:12.90 


+31:23:29.5 


LAL 258 


II 0.07 


-1.52 


105 


03:29:13.04 


+31:17:38.4 


ASR 28; LAL 263 


II 0.00 


-1.17 



Table 3 — Continued 



Index 




jjeci. 


Lit. Matches 




A c 
Ak 


aiRAC 


106 


03:29:13.13 


+31:22:52.9 


LAL 262 


jj 


1.02 


-1.30 


107 


03:29:15.34 


+31:29:34.6 




jj 




-0.39 


108 


03:29:15.66 


+31:19:11.1 


ASR 59 


JJ 




-1.08 


109 


03:29:16.59 


+31:23:49.6 


LAL 276 


JJ 


0.81 


-0.70 


110 


03:29:16.81 


+31:23:25.4 


LAL 279 


JJ 


0.24 


-0.48 


111 


03:29:17.66 


+31:22:45.2 


SVS 2; LAL 283 


JJ 


0.31 


-0.81 


112 


03:29:17.76 


+31:19:48.2 


ASR 80; LAL 286 


JJ 


0.40 


-1.33 


113 


03:29:18.65 


+31:20:17.9 


ASR 81 


JJ 




0.00 


114 


03:29:18.72 


+31:23:25.5 


LAL 293 


JJ 


0.08 


0.08 


115 


03:29:20.05 


+31:24:07.6 


LAL 296 


JJ 




0.13 


116 


03:29:20.42 


+31:18:34.3 


SVS 5; IRAS 03262+3108; ASR 112; LAL 300 


JJ 


1.55 


-0.60 


117 


03:29:21.30 


+31:23:46.5 


LAL 302 


JJ 




-1.02 


118 


03:29:21.55 


+31:21:10.5 


LAL 304 


JJ 


0.00 


-1.25 


119 


03:29:21.87 


+31:15:36.4 


ASR 123; LAL 307 


JJ 


0.21 


-1.49 


120 


03:29:23.15 


+31:20:30.5 


ASR 78; LAL 310 


JJ 


0.00 


-1.02 


121 


03:29:23.23 


+31:26:53.1 


LAL 308 


JJ 


0.16 


-1.02 


122 


03:29:24.08 


+31:19:57.8 


ASR 79; LAL 313 


JJ 


0.58 


-0.39 


123 


03:29:25.92 


+31:26:40.1 


SVS 4; LAL 318 


JJ 




0.06 


124 


03:29:28.01 


+31:25:11.0 




JJ 




-1.00 


125 


03:29:29.79 


+31:21:02.8 


LAL 333 


JJ 


0.39 


-0.81 


126 


03:29:30.39 


+31:19:03.4 


LAL 336 


II 


0.00 


-1.05 


127 


03:29:32.56 


+31:24:37.0 


LAL 342 




0.38 


-1.43 


128 


03:29:32.86 


+31:27:12.7 


LAL 344 




0.08 


-1.67 


129 


03:29:35.70 


+31:21:08.6 








-0.27 


130 


03:29:37.63 


+31:02:49.3 






0.00 


-1.14 


131 


03:29:37.72 


+31:22:02.6 






0.00 


-1.29 


132 


03:29:44.15 


+31:19:47.5 








-0.64 


133 


03:29:54.03 


+31:20:53.1 






0.00 


-0.99 


134 


03:29:02.69 


+31:19:05.8 


ASR 62; LAL 187 


II/III 




-1.71 


135 


03:29:09.41 


+31:14:14.1 


ASR 54 


II/III 




-0.32 


136 


03:29:26.79 


+31:26:47.7 


SVS 6; LAL 321 


II/III 


0.22 


-2.14 


137 


03:29:29.26 


+31:18:34.8 


LAL 331 


II/III 


0.65 


-2.12 



''J2000 

''Asterisks mark "deeply embedded" sources with questionable IRAC colors or incomplete IRAC photometry and 



relatively bright MIPS 24 nm photometry. 

"^Only provided for sources with valid JHKg photometry. 

"^Extinction is not accounted for in these values. High extinction can bias ajuAC to higher values. 
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Table 4. Summary of Disk Fraction Measurements by Region 



Region 


R.A.» 


Decl.^- 


Radius (arcsec) 


No. with IR-excess 


Total No. of Members 


Disk Fraction 


Main 


03:29:04.67 


+31:19:57.9 


330 


72 


87 


0.83 + 0.11 


SI 


03:28:55.12 


+31:16:23.8 


90 


11 


15 


0.73 ± 0.07 


S2 


03:29:03.47 


+31:16:54.4 


90 


6 


7 


0.86 + 0.10 


Nl 


03:29:10.64 


+31:23:01.4 


90 


13 


14 


0.93 + 0.07 


N2 


03:29:02.28 


+31:21:29.7 


90 


14 


16 


0.88 + 0.07 


E 


03:29:27.34 


+31:19:57.7 


120 


7 


7 


1.00 + 0.13 



''J2000 
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Fig. 1. — Color-composite image of the IRAC mosaics of NGC 1333; 3.6, 4.5, and 8.0 /im 
images are mapped to blue, green, and red, respectively. 
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Fig. 2. — Color-composite image of the IRAC and MIPS mosaics of NGC 1333; 4.5, 8.0, and 
24 /im images are mapped to blue, green, and red, respectively. 




Fig. 3. — Color-color diagrams used for identifying PAH emission sources. The upper plot 
is the [4.5] — [5.8] vs. [5.8] — [8.0] color-color diagram, and the lower one is the [4.5] — [8.0] 
vs. [3.6] — [5.8] color-color diagram. PAH-emission sources are marked with crosses, and 
broad-line AGN are marked with plusses. 




Fig. 4. — [4.5] vs. [4.5] — [8.0] color-magnitude diagram, used for identifying likely broad-line 
AGN. PAH-emission sources are marked with crosses, and broad-line AGN are marked with 
plusses. 
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Fig. 5. — Spatial distribution of IR excess sources in the NGC 1333 region. At left are those 
sources flagged as hkely extragalactic contaminants. Crosses mark PAH-emission sources 
and plusses mark likely broad-line AGN. Note that their distribution is relatively uniform 
and has httle to do with the cluster at the center of the field of view, nor the extinction 
distribution overlaid in gray contours {Ay — 3, 5, 7, 9, 11, 13, 15). At right are the identified 
YSO sources, overplotted on the same extinction contours. Crosses mark all Class I sources 
and plusses mark the more-evolved Class II sources. The clustered distribution of these 
sources and their close association with the extinction distribution are distinct. 
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Fig. 6. — [5.8] magnitude histogram comparing numbers of sources filtered as contaminants 
compared to the model prediction derived from the Bootes field analysis presented in Ap- 
pendix \M The top histogram is for the PAH emission sources, and the bottom is for the 
likely broad-line AGN. Bootes field models, scaled for the NGC 1333 coverage area and ad- 
justed for relative areas of differing extinction, are plotted as black dot-dashed lines. The 
distribution of sources filtered from the NGC 1333 source list are plotted as dashed lines. 
The distribution of all IR excess sources in NGC 1333, excluding the PAH emission sources, 
is plotted as the gray solid line. The parallel slanted line shading marks the region of the 
plot for sources dimmer than the 90% completeness limit of [5.8] = 14.5. 
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Fig. 7. — Color-color diagrams used for YSO identification and classification. The left plot 
is the [3.6] — [4.5] vs. [4.5] — [5.8] color-color diagram, used for identifying likely proto- 
stars and filtering out barely-resolved shock emission. Diamonds mark the IRAC-identified 
protostellar sources, and circles mark the more-evolved Class II sources. Triangles mark 
transition/debris disk candidates, sources with colors consistent with reddened photospheres 
that have IR-excess at 24 //m. Asterisks mark hkely shock emission knots. Squares mark 
"deeply embedded YSOs" , hkely protostars, that are composed of shock knot and broad-line 
AGN candidates and sources with incomplete IRAC four-band photometry that have partic- 
ularly red IRAC/MIPS colors. The right plot is the [K — [3.6]]o vs. [[3.6] — [4.5]]o color-color 
diagram, used for secondary YSO identification and classification for sources that lack either 
5.8 or 8.0 ixm photometry. 
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Fig. 8. — Two IRAC four-band color-color diagrams used for classifying YSOs. The top plot 
is the [3.6] — [4.5] vs . [5.8] — [8.0] colo r-color diagram, used in several previous papers (e.g. 
Megeath et Zlbooi lAUen et all 120041 ). The bottom plot is the [3.6] - [5.8] vs. [4.5] - [8.0] 
color-color diagram, used in this work for identifying likely Class II sources. Plot overlays 
are the same as Figure [71 
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Fig. 9. — [4.5] — [5.8] vs. [5.8] — [24] color-color diagram, used for identifying likely transi- 
tion/debris disk sources and to verify protostellar classifications. Plot overlays are the same 
as Figure [3 
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Fig. 10. — SCUBA 850 fim contours with Class II sources overlaid on the left in green and 
Class I sources on the right in red. The first contour level marks 0.05 Jy/beam, and each 
successive contour is double the previous level. Note that the Class I population traces the 
dense material much more closely than the Class II sources. 
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Fig. 11. — Normalized histograms comparing the SCUBA 850 fim flux at the positions of 
Class I and Class II populations in gray and black, respectively. This flgure demonstrates 
that protostars are much more strongly conflned to regions of high gas density compared to 
their more evolved counterparts. 
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Fig. 12. — Histograms of nearest neighbor distances in NGC 1333. Top and middle his- 
tograms are the distributions among the two evolutionary classes, protostars and Class II 
respectively. The bottom histogram is the distribTition for all YSOs together. The dashed 
vertical lines mark the range where all three distributions have clear peaks. 
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Fig. 13. — A comparison of NNn surface density mapping results for the NGC 1333 YSOs, 
using n = 4, 6, 11, and 18 neighbors. As the number of neighbors used is increased, smaller 
scale structure is overwhelmed by the large scale distribution that is dominated by the 
"double c luster". Contours leve ls mark stellar surface densities of 1 cr below successive 
contours ( ICasertano fc Hutlll985l ). and the linear inverse grayscale sets 1000 pc~^ as black 
and pc~^ as white. 
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Fig. 14. — At left: NN6 map of the surface density of all sources detected at Ks, [3.6] 
and [4.5] /xm with Ks < 15. At right: NN6 map of the surface density of 5]5iteer-identified 
YSOs only, i.e. only confirmed cluster members. The circle labeled "Main" marks the cluster 
region used for measuring the fraction of members with circumstellar material. Other smaller 
regions are also examined and their results are listed in Table HI The box marks the control 
field used for characterizing field star contamination in the statistically inferred Class III 
diskless pre-main sequence star population. 
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Fig. 15. — Ks magnitude histogram (KMH) demonstrating the measurement of the total 
membership of the cluster region in Figure [TH The top plot shows the measured KMH for 
those sources within the circle that are detected at H, Ks, [3.6] and [4.5] /im with Ks < 14. 
The middle plot shows the KMH of the control field, scaled by the ratio of the areas of the 
cluster region and the control field and shifted to account for the mean extinction difference 
between the two fields. The bottom plot shows the final field star corrected cluster KMH 
for the cluster region, resulting in a total membership estimate of 82. Overplotted in red is 
the KMH for all S^^rfzer-identified YSOs detected at Ks band in this region. There are 72 
with Ks < 14, thus the fraction of members with disks is 83%. 
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Fig. 16. — Stellar surface density as a function of radius for the 5'pifeer-identified YSOs 
centered on four distinct choices of cluster center points. In the left column of plots, the 
straight, black, dashed lines are power laws of a = (shallow slope) and a — —3.3 (steep 
slope) that were chosen to fit the functional behavior of the radial profiles at small and large 
radii, respectively. Also in the left column, a function of the form a — (7o(l + {—YY^ 
plotted in solid gray, with parameters ctq = 220 pc^^, /9 = 6, 7 = 3.3, and ro = 0.3 pc, the 
location of the knee in the two power law model. In the right column, the overlays are the 
best fit King profile in dashed black and the best fit EFF profile in solid gray for the "Main" 
radial profile. 
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Fig. 17. — YSO number counts per 45° azimuthal bin (Nyquist sampled), with error bars 
derived from Poisson counting statistics for the contents of each bin. The dashed hne 
marks the mean source count per bin. The AAP is the standard deviation of the values 
for each bin about this mean, divided by the square root of the mean. The resulting AAP 
of 2 for NGC 1333 is quite high in reference to other young clusters analyzed this way 
(iGutermuth et al.ll2005l ). The substructure that is dominating this measurement is likely 
the two-peaked nature of the "double cluster" . 
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Fig. 18. — Color-color diagrams of the Bootes field data used to identify PAH emission 
sources. The top plot is the [4.5] — [5.8] vs. [5.8] — [8.0] color-color diagram, and the botton 
one is the [3.6] — [5.8] vs. [4.5] — [8.0] diagram. Diamonds mark residual contaminants 
that have mid-IR colors and magnitudes similar to Class I YSOs. Circles mark residual 
contaminants with Class II colors and magnitudes. 




Fig. 19. — [4.5] vs. [4.5] — [8.0] color-magnitude diagram, used for identifying likely broad-line 
AGN. Diamonds mark residual contaminants that have mid-IR colors and magnitudes similar 
to Class I YSOs. Circles mark residual contaminants with Class II colors and magnitudes. 
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Fig. 20. — Source counts per square degree vs. [5.8] for the Bootes field. The upper left 
histogram is the distribution of PAH emission sources. The upper right solid black histogram 
is all other sources with IR-excess emission in the Bootes field, assumed to be dominated 
by broad-line AGN. The gray dot-dashed histrogram is the distribution of sources that are 
successfully filtered with our broad-line AGN cut shown in Figure [191 Gray solid straight 
lines in both plots are the power law fits presented in the text. The residual contaminating 
sources are classed as they would be if they were YSOs to gauge contamination by class. 
The lower left histrogram is for residual "protostar" contaminants, while the lower right is 
for residual "Class 11" contaminants. 



